RADIATIVE DECAY OF BUBBLE OSCILLATIONS IN A 
COMPRESSIBLE FLUID 

A. M. SHAPIRO* AND M. I. WEINSTEINt 

Abstract. Consider the dynamics of a gas bubble in an inviscid, compressible liquid with 
surface tension. Kinematic and dynamic boundary conditions couple the bubble surface deformation 
dynamics with the dynamics of waves in the fluid. This system has a spherical equilibrium state, 
resulting from the balance of the pressure at infinity and the gas pressure within the bubble. We study 
the linearized dynamics about this equilibrium state in a center of mass frame. We prove that the 
velocity potential and bubble surface perturbation satisfy point- wise in space exponential time-decay 
estimates. The time-decay rate is governed by the imaginary parts scattering resonances. These are 
characterized by a non-selfadjoint spectral problem or as pole singularities in the lower half plane of 
the analytic continuation of a resolvent operator from the upper half plane, across the real axis into 
the lower half plane. The time-decay estimates are a consequence of resonance mode expansions for 
the velocity potential and bubble surface perturbations. The weakly compressible case (small Mach 
number, e), is a singular perturbation of the incompressible limit. The scattering resonances which 
govern the remarkably slow time-decay, are Rayleigh resonances, associated with capillary waves, 
due to surface tension, on the bubble surface, which impart their energy slowly to the unbounded 
fluid. Rigorous results, asymptotics and high-precision numerical studies, indicate that the Rayleigh 
resonances which are closest to the real axis satisfy = C'(exp(— k We , k > 0. Here, 

We denotes the Weber number, a dimcnsionlcss ratio comparing inertia and surface tension. To 
obtain the above results we prove a general result estimating the Neumann to Dirichlet map for the 
wave equation, exterior to a sphere. 

Key words. 

AMS subject classifications. 

1. Introduction. Consider a gas bubble, surrounded by an unbounded, inviscid 
and incompressible fluid with surface tension. This system has a family of equilibria, 
consisting of translates of a spherical gas bubble, whose radius is set by a balance of 
pressure at infinity with pressure inside the bubble. In this paper we prove pointwise 
in space time-decay estimates for the linearized evolution near this family of equilibria. 
We also obtain very precise information on the rate of decay. 

Background: The dynamics of gas bubbles in a liquid play an important role in many 
fields. Examples include underwater explosion bubbles [10| (15 cm), seismic wave- 
producing bubbles in magma [30) . bubbly flows behind ships and propellors (1 cm), 
bubbles at the ocean surface [H] (0.015 — 0.5 cm), microfluidics [46] (50 /im), bubbles 
used as contrast agents in medical imaging [9] (2 ^m), and sonolumincscence [3|l27j 
(0.1 — 10 ^m). For a discussion of these and other applications of bubble dynamics, 
see the excellent review articles [TnHH] and references cited therein. 

The dynamics of a bubble in a liquid are governed by the compressible Navier- 
Stokes equations in the liquid external to the bubble, a description of the gas within 
the bubble, and boundary conditions (kinematic and dynamic) which couple the fluid 
and gas. We assume the gas inside the bubble to be at a uniform pressure throughout 
and to satisfy a thermodynamic law relating the bubble pressure to the bubble volume. 
Rayleigh initiated the study of bubble dynamics and derived an equation for the radial 
oscillations of a spherically symmetric gas bubble in an incompressible, inviscid liquid 
with surface tension |131I28] . This problem has a spherical equilibrium, balancing the 
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pressure at infinity and the pressure within the bubble. A general (asymmetric) per- 
turbation of this spherical equilibrium will excite all harmonics and Rayleigh showed, 
in the linear approximation, that each harmonic executes undamped time-periodic os- 
cillations |13| . With viscosity present, the pulsating bubble oscillations are damped 
and the spherical equilibrium is approached. 

A second very important damping mechanism, energy preserving in nature, is due 
to compressibility of the fluid. In a non-viscous and compressible fluid, a perturbation 
of the equilibrium bubble will, due to coupling at the gas-liquid interface, generate 
acoustic waves in the fluid which, in an unbounded region, propagate to infinity. The 
bubble dynamics acquire an effective damping, due to energy transfer to the fiuid, and 
its propagation to infinity. If surface tension is included asymmetric modes should 
damp and the bubble shape should approach that of a sphere as time advances. 

Keller and co-workers [Sl IlOlfTT] modeled slight compressibility of the fiuid by the 
linear wave equation, exterior to the bubble. Both the bubble interior pressure law 
and the boundary conditions are kept as in the incompressible case. In the spherically 
symmetric setting this leads to an ODE, which captures acoustic radiation damping. 
A systematic derivation of the model of Keller et al., in the spherically symmetric 
setting, was presented by Prosperretti-Lezzi [TTll^ . 

We expect, for the inviscid, compressible problem with surface tension, that the 
family of translates of the spherical equilibrium is (locally) nonlinearly asymptotically 
stable. Specifically, a small perturbation of the equilibrium spherical bubble, will 
induce translational motion of the bubble and deformation of its surface. We expect 
that, in a frame of reference which moves with the bubble center of mass, ^cm(Oj ^ 
small perturbation of the spherical bubble will damp toward a spherical equilibrium 
shape. As in many studies of asymptotic stability for coherent structures in spatially- 
extended conservative nonlinear PDEs, linear decay estimates of the type established 
in this work can be expected to play a role in the estimation of the convergence to 
equilibrium for the nonlinear dynamics. 

1.1. PDEs for a gas bubble in a compressible liquid. Consider a gas bub- 
ble, occupying a bounded region B{t), surrounded by an inviscid and compressible 
fiuid with surface tension. We consider the boundary of the bubble, dB{t), to be 
parametrized by a function R : {a,t) i-> R(a,t) G dB{t) C R'^, with parameter a, 
e.g. spherical coordinates. 

The time-evolution of the fluid is governed by the system: 

5ju+ (u- V)u+ ivp(p) = 0, xeM'''\B(t) (1.1a) 
P 

dtp+V ■ (pu) =0, x£R^\ B{t) (1.1b) 

(uo R) • n = 9tR • n, dB{t) (1.1c) 

Pbubbiolas(i) - Pfiuid|aB(4) = 2tTiJ[R], dB{t), (l.ld) 

where H[R] = i V • fi denotes the mean curvature at location R G dB{t)-. see, for 
example, [T3| (Article 275, Equation 5), or [7] (Section 3-3). 

The pressure within the fiuid is assumed to obey an equation of state: p = 
Pfluid = p{p)- Equations (jl.lap and (jl.lb[) express conservation of momentum and 
mass. Equation (jl.lcp is the kinematic boundary condition, i.e. the normal velocity, 
SfR • n, of the material point on the bubble surface moves with the normal velocity of 
the fiuid. The Young-Laplace boundary condition, also called the dynamic boundary 
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condition, (jl.ld[) . expresses that the jump in pressure at the fluid bubble interface is 
proportional to the mean curvature |15] : 

Phuhhic\gB(t) ^ MuidlaB(t) ^ Surface Tension x (2 x Mean Curvature). (1.2) 

We assume that the pressure within the bubble is spatially uniform and given by the 
polytropic gas law [19]: 



PbubbIoUw=^B = ^g(^=Pc,(^P(^j , 7>1- (1-3) 

The pressure within the fluid is assumed to satisfy a general smooth relation of the 
form: p = p{p), so that the system is determined by only one state variable. 

Finally we assume that the initial velocity is irrotational, V A Uq = 0. It fol- 
lows that the velocity field remains irrotational for all t > 0; see, for example, [T^ . 
Article 33. Thus, there is a single- valued velocity potential $, such that u(x, t) = V$. 

Equilibrium solutions: Equations p.ip have a spherically symmetric equilibrium 
solution: 

u = 0, R = ar= a- p = Poc, P^Poo- (1-4) 

|x - Xo| 

The equilibrium bubble radius, a, is uniquely determined via the dynamic boundary 
condition 

k 2(7 2(7 , , 

Ti ^ Poo = — or Pea = Poo H • (1-5) 

(^a^y a cq /^oo ^ y , 

We consider these dynamics in the linear approximation. Introduce spherical 
coordinates, (r, 9, (f) = (r, 51), with the origin chosen to be the bubble center of mass, 
^cm(^)- We express a small perturbation of the spherical bubble as: 

$(r,ri,0 = constant + ^(r,^!,^, r=|x-|,^(0|, (1.6) 

In Appendix [Clwe re- write the system (|l.ip . relative to coordinates centered at Ccm(^)- 
In the linear approximation, ^crn(^) ^ ^cmlofl and the nondimensional system of 
equations, linearized about the spherical equilibrium bubble are: 

e^^t^* - A* = r > 1, (1.8a) 

*r - A r = 1, 

(1.8b) 

(/3,rr)i.(5.) = o. H<i 

(1.8d) 



^At nonlinear order, ^(.m(*) typically evolve. 
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Here, e denotes the Mach number (m = e is used in the derivation of the non- 
dimensional equations in Appendix [C|), Ca, the Cavitation number, and We, the 
Weber number. 

We shall focus on the initial-boundary value problem for (|1.8p with data corre- 
sponding to an initial perturbation of only the bubble surface: 

^{r, n,t = 0) = dt^{r, 17, t = 0) = 0, /3(t = 0, O) given and sufficiently smooth. 

(1.9) 

A5 denotes the Laplacian on S*^, in spherical coordinates given by (|A.6p . The spherical 
harmonics, Y[^{yi), are eigenfunctions: — AsYJ™ = l{l~\-\)Yj^ , I >Q, < Z, forming 
a complete orthonormal set in L'^{S^) with respect to the inner product (ai, 02) l^(s'^)', 
see Section [Dl The orthogonality conditions (jl.8d[) . derived in Appendix [UJ express 
our choice of coordinates (in the linearized approximation) placing the origin at the 
bubble center of mass. 

1.2. Overview of results and discussion. We conclude this section with an 
overview and discussion of results. 

1. Time-decay of solutions to wave equation on M''\{|x| < 1} with time- 
dependent Neumann data: Theorcm l4.2l is a general result, of independent interest, 
on the time-decay and resonance expansion of solutions to the initial-boundary value 
problem for the wave equation on — 5^. The data prescribed on are assumed 
sufRciently smooth and exponentially decaying with time. Theorem 14.21 generalizes 
the results of glllll]; see also [Tl|39l. 

2. Exponential time-decay estimates for the bubble surface perturba- 
tion (Theorem lS.ip : Theorcm l4.2l on the Neumann to Dirichlet map, together with 
the detailed information we obtain on the locations of scattering resonances in the 
lower half plane, is used to prove that the solution of the initial value problem (|1.8p 
with initial data (|1.9p decays exponentially to zero, pointwise, at a rate C'(e~'"^'^^*), 
r(e) > 0, as t tends to infinity. Moreover, the linearized velocity potential, ^{r,^l,t) 
and bubble surface perturbation, /3{^l, t), satisfy resonance expansions in terms of out- 
wardly radiating states of the scattering resonance problem. The expansion converges 
in C^{K X M+), where K denotes any compact subset of |x| > 1. 

3. Scattering resonances and radiation damping: The exponential rate 
of decay, r(e) = |3Aj(e)|, is determined via the scattering resonance problem^ (j3.4p . 
a non-selfadjoint spectral problem associated with the time-harmonic solutions of the 
linearized compressible Euler equations -I- boundary conditions on |x| = 1. with outgo- 
ing radiation conditions at infinity. Two families of scattering resonances, associated 
with the Helmholtz equation A^E* + w^^P ~ for |x| > 1, appear in the linearized 
problem. Rigid resonances, {w(,fc(e)};>o,i<fe<(+i fThcorcm 13. ip . are associated with 
Neumann [sound soft) boundary conditions imposed on the sphere, and outgoing ra- 
diation conditions at infinity and deformation resonances, {A/,fc(e)}/>o,i<fc<z-i-2 (The- 
orem 16. ip . associated with the hydrodynamic boundary conditions on the sphere, 
responsible for the deformation of the bubble-fluid interface, and outgoing radiation 
conditions at infinity. Theorem 16.31 implies that for each fixed e > 0, there is a 
strip containing the real axis, in which there are no scattering resonances; hence 
|9A±(e)| >0. 

4. Asymptotics of deformation resonances and Fermi's Golden Rule: 

In the incompressible limit, e 4, 0, all resonances (rigid and deformation) have imag- 
inary parts which tend to minus infinity, except for the sub-family of deformation 
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resonances ealled Rayleigh resonances, {Aj*(e)}/>o- Raylcigh resonanees are scatter- 
ing resonances in the lower half plane to which the real Rayleigh eigenfrequencies (of 
undamped oscillations in the incompressible problem, e = 0) perturb upon inclusion 
of small compressibility, e > 0. Their detailed asymptotics for e small is given in 
Theorem 16.21 In particular, we find for the 
Imaginary parts of the Rayleigh Deformation Resonances: 



i(^Ca + 2(3.-l)4 



i[(/ + 2)(.-l)r(. + iy[|i] ( 



Wei ^^'ivWc, 



(1.10) 



for I = 2,3,... . The proof of (jl.lOp relies on use of detailed properties of spherical 
Hankel function and, in particular, a subtle result on the Taylor expansion of the func- 



tion z i-> Gi{z) = zdh\^\z)/h\'->{z) in a neighborhood of z = 0; sec ProDOsition l8.4l 

The infinite set of pairs of (real) Rayleigh eigenfrequencies {A^ (0)}i=o,2.... may be 
viewed as: embedded eigenvalues in the continuous spectrum of the unperturbed (e = 0, 
incompressible) spectral problem. The negative imaginary part in (|1.10p is an instance 
of the Fermi Golden Rule. An expression coined originally in the context quantum 
electrodynamics [11 , it refers to the induced damping of an "excited state" (the bubble 
perturbation), due to coupling of an "atom" (the deforming bubble) to a field (wave 
equation); see also, for example, [29 l l33 l l44] . 

5. Scaling behavior of the decay rate, |5Aj(e)| (Section [7]): Theorem l6.21 
asymptotics and high-precision numerics show that the exponential rate of decay is 
given by |5A^(e)| — |3A~(e)| > 0, where Aj(e) are the scattering resonance energies 
closest to the real axis, and in the lower half plane. These resonances Aj(e) ~ 
± ; _ (e'^We) . Moreover, an asymptotic study of the results of 



,(1 



A^(^-|, where h{e) 
Theorem 16.21 yields: 



3?a: 



^^A±(6) 



■O 



/We 



n>0. 



In contrast, the monopole (spherically symmetric) resonance has an imaginary part 
which isO{e). This remarkably slow rate of decay is related to the scattering resonance 
problem being a singular perturbation problem in the (incompressible) limit of small 
e; the wave equation e^df^ = A$ reduces to Laplace's equation A$ = as e — > 0. 
Figure 11.11 displays, for a particular small choice of e, a range of resonance energies 
including those located closest to the real axis. Physically, these very slowly decaying 
bubble shape modes are associated with capillary waves on the bubble surface, are 
excited only by asymmetric perturbations, and very slowly transfer their energy to 
the infinite fiuid. □ 



■^We remark on an interesting example, where scattering resonances converge to the real axis. 
Stefanov and Vodev 1361 show, for the system of linear elasticity on the exterior of a ball in K.^, that 
the scattering resonances converge to the real axis exponentially fast, as the spherical harmonic index, 
I, tends to infinity. Consequently, the time-decay of solutions is not exponentially fast; see [4l l35ll38] 
and references cited therein. The modes associated with the resonances which come ever closer to 
the real axis are called Rayleigh surface waves, a decay rate limiting mechanism analogous to the 
capillary surface waves {Rayleigh resonance modes) of the bubble problem wc consider in this article. 
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Fig. 1.1. (a) Numerically computed scattering resonances associated with spherical har- 
monics I = 20, 24, 28, . . . , 52 for e = 0.1, We = 1. Rayleigh resonances are those located just 
below the real axis, (h) Vertical axis is scaled to show the resonances closest to the real axis. 



6. Monopole vs. Multipole radiation: A discussion of emission of acous- 
tic radiation in the physics literature (see jl9| and references cited therein) is based 
on a heuristic energy balance argument, which assumes that the dominant emis- 
sion of acoustic radiation is through the monopole {I = 0, purely radial) mode. 
This assumption implies decay of the bubble perturbation on the time scale of or- 
der Tnionopoic(e) ~' e^^. Our results demonstrate that there arc indeed non-symmetric 
vibrational modes of a bubble, which have the much longer lifetime 

T(e) - e (e"^Wc) cxp {k f-^Wc) > T,„onopoio(e), k > 0. 

Of course, in many systems one would expect viscosity (not included in the modeling) 
to contribute a dominant correction to the imaginary part. 

We note a relation of these results to those in [18] and [21], who used energy flux 
arguments to study the radiation of acoustic energy in compressible atmospheres. 

Some future directions: In addition to the question of nonlinear asymptotic sta- 
bility of the family of spherical equilibrium bubbles, we remark that there is a rich 
class of solutions in the radially symmetric setting. In the case of radial symmetry, 
Rayleigh showed that the dynamics of a spherical bubble in an incompressible liquid 
exactly reduce to a nonlinear ordinary differential equation for the time-evolution of 
the radius. This equation admits the spherical equilibrium (constant radius) state as 
well as time-periodic (radially pulsating) states. Plesset [21] extended the analysis 
in the spherically symmetric case to include viscosity. It is of interest to study the 
dynamics near the spherically symmetric time-periodic states of the Rayleigh-Plesset 
equations. We expect that these would be unstable; periodic oscillations would couple 
to continuous spectral modes, resulting in radiation-damped "breather" oscillations; 
see, for example, |34j and references cited therein. The period map or monodromy 
operator associated with the linearization about such a state, for e = 0, would have an 
embedded Floquet multiplier on the unit circle, which would perturb to an unstable 
Floquet multiplier outside the unit circle [201132] . 

Finally, we remark that the problem we consider is one of a very large class, 
involving an infinite dimensional conservative system comprised of two coupled sub- 
systems: one subsystem has a discrete number of degrees of freedom (here the mode 
amplitudes of the spherical harmonics of the bubble surface) and one has a continuum 
of degrees of freedom, the velocity potential governed by the wave equation. While 



RADIATIVE DECAY OF BUBBLE OSCILLATIONS IN A COMPRESSIBLE FLUID 



7 



there has been significant progress on such systems, where the number of "sohton" 
degrees of freedom of the discrete subsystem is finite (for example, see [T^[M1H5] and 
references cited therein), systems hke the bubble-fluid system, which involve the cou- 
pling of infinite dimensional systems, are central in physical and engineering science 
and are an important future direction. 

Acknowledgements: The authors thank D. Attinger, J.B. Keller, A. Soffer, E. 
Spiegel, J. Xu and M. Zworski for stimulating discussions. This research was sup- 
ported in part by US National Science Foundation grant DMS-07-07850 and NSF 
grant DGE-0221041. The Ph.D. thesis research of the first author was supported, in 
part, by the IGERT Joint Program in Applied Mathematics and Earth and Environ- 
mental Science at Columbia University. Some of the discussion of this paper, appears 
in greater detail in [31j . 

1.3. Some oft used notation. 

• (a, /?) = (a, /3)i2(g2') = inner product on L^(iS'^); Section [D] 

• As = Laplacian on 5"^; Equation (|A.6|) . Yj^iO, (f)) = ^"'(ri), spherical har- 
monic; Appendix [Bl 

• Dimensionless quantities: We = Weber number, Ca — Cavitation number, 
e = M = Mach number; Appendix [Cl 

,,[fCa + 2(3,-1)^, /=.0 ^^^^^^ 



We 



{l + 2){l-l), l>l, 



2. Conservation of Energy. An important role in our analysis of time-decay 
of solutions of the initial-boundary value problem for (|1.8p is played by the following 
conservation law. We use it to establish the scattering resonance frequencies as lying 
in the lower half plane; see the proof of Proposition 16. II 

Proposition 2.1. Let ( 4'(x,t), /3(ri,i) ) denote a smooth complex-valued solu- 
tion of the linearized equations. 

(i) Then if ^ decays sufficiently rapidly as \x\ — > oo, then the following func- 
tional is time-invariant 

£ = f (e2|a,*|2 + |vvl/|2)dx+37(f + 4)|(/3'^o")|'+w?((-As-2)/3,/3), (2.1) 

J\x\>l 

where {oL,j3) denotes the inner product on L'^{S'^). 

(ii) The energy, £, can also be expressed as: 

£= [ {e^\d,^\^ + \V^\')dx + (fCa-f 2(37-1)^) \{f3X)\' 

J\x\>l 



4 E E(^ + 2)^-1) iw^ynf- (2.2) 

l>2 \m\<l 



Furthermore, £ is positive definite on solutions of the linearized equations (|1.8p since 
7 > 1. 

Proof. Part (i): Multiplying the wave equation (|1.8al) by dt'^ and taking the real 
part of the resulting equation yields 

9te(^',at*) + v-(-a^v* - at*W) = 0, (2.3) 
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where e(5',94^') — e^|9t4'p + |V\l/p. Integration of the local conservation law 
over the region > 1 yields 

^ / e{^,dt-^) dx+ [ {d^V-^ ■n+ dt'HV^ -h) dS ^Q, (2.4) 

where n denotes the unit normal to 5^, pointing into the ball. The theorem no follows 
from re-expressing the second term in (|2.4p using the boundary conditions: 

(9^V* •n + af^'W-n) dS 



\x\ = l 

( ^ -4- ^ 1 I (/-(. 1^ 

dt 



37(^ + 4)|(/3,roO)|^ + ^((-As-2)/?,/3) (2.5) 



Part (ii): We show now, using the center of mass constraint, (jl.Sdp . that the energy 
is positive definite on solutions of the linearized system, (|1.8p . 

By (jl.8d[) we can expand /? in spherical harmonics as /? = {P^Y^) + /3;>2, 
where pi>2 = Ei>2E|™|</ (/^-^D ^f"- Consider the expression ((-A5 - 2)/?,^). 
We have, by orthogonality 

((-As - 2)/3,/3) = -2 + ((-A5 - 2)A>2, /3/>2) 

= -2 |(/?,ro°)r + E E + 2)(^ - 1) l(A>^r>l' (2.6) 

l>2 \m\<l 

Use of (|2.6p in the expression for the energy (|2.ip yields (j2.2p , which is clearly positive 
definite. □ 

3. The Scattering Resonance Problem. Our results on time-decay of the 
velocity potential and bubble surface perturbations (Theorems 14.21 and 15 . 1 p are proved 
using their representation as inverse Laplace transforms of a resolvent operator applied 
to the initial data; see Equations (j4.19p and (|5.12p . in which the contour of integration 
is a horizontal line in the upper half of the complex frequency plane. Time decay 
(and a resonance mode expansion) is deduced by analytically continuing the resolvent 
kernel (Green's function) and deforming the integration contour across the real axis 
(essential spectrum) into the lower half plane. The time-decay is determined by 
poles of this analytically continued resolvent kernel. These poles are called scattering 
resonances or scattering poles. An alternative characterization of scattering resonances 
is as complex frequencies in the lower half plane for which there are non-trivial time- 
harmonic solutions of the linearized equations, satisfying an outgoing radiation (non- 
self-adjoint) boundary condition at infinity. 

In this section we use the alternative (spectral) characterization as solutions of a 
non-sclfadjoint eigenvalue problem. We consider two classes of scattering resonances: 

1. Rigid resonances, associated with the wave equation on the exterior of the 
(rigid) unit sphere with Neumann boundary conditions, and 

2. Deformation resonances, associated with the linearized system (|1.8|) . 

Both families of resonances, for e > 0, lie in the lower half plane. Of particular 
interest is the resonance in each family with the smallest (in magnitude) imaginary 
part. Denote by w*(e) and A^(ejl the rigid and deformation resonances of smallest 
imaginary parts. We find for e small: 

|3c^4e)| = O (e-i) , \'^K{e)\ = ©(f-^e"*), k>Q ^ I^A^e)] « l^cj^e)] . (3.1) 



^Here, and in later sections, A^, (e) may be denoted simply A*(e). 
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3.1. Rigid - Neumann scattering resonances. We consider the scattering 
resonance problem, associated with the wave equation e^df"^ = A^* in the exterior of 
the unit sphere with Neumann boundary conditions. Seeking time-harmonic solutions, 
^'(7-, f2, t) = e~*"**ij, which are outgoing as r — )■ 00 , we arrive the following eigenvalue 
problem 

( A + (ew)2 )^'^(r,l)) =0, r>l; 9,*^(r, 17) - 0, r = 1. (3.2) 

Outgoing solutions are spanned by functions of the form h[^\eu!r) yj'"(51). Thus, u) 
is a Neumann scattering resonance if and only if dhi{euj) = . The following result 
summarizes results in [22l|4T]; see also [1]. 

Theorem 3.1. Fix an arbitrary e > and arbitrary. 

(i) For each I > 0, the equation (ew) = has a family of solutions 
{w;,fc(e) = e^^w/^fc : fc = 1, . . . , Z + 1}. 

(ii) The set of all Neumann scattering resonances is a discrete subset of the lower 
half complex plane and is uniformly bounded away from the real axis. In particular, 
there exists > 0, 1 < ^* < ^* + 1 such that a;*(e) =dcf (e) is a resonance whose 
imaginary part is of minimal magnitude, i.e. 

9w,,fc(e) < 3?a;4e) < 0, for all I >0, k = l,...,l + l (3.3) 

(iii) There exist constants Ci, C2 > 0, such that for all I > 

Ci e-Mi/3 < |3^,,^(e)|, /»i and |c^/,fe(e)| = e^^M < C2e"M . 

This result is essentially due to Tokita (41] and uses fundamental results on the 
asymptotics of Hankel and Airy functions; see Olver [Tll22|. Explicit approximations 
to zeros of dH^ are given in Equation (|E.17p . 

3.2. Deformation resonances. We seek time-harmonic solutions of the lin- 
earized perturbation equations ([L5)) : * = e"'^*\E'A(r, O), (3 = e^*^*/3xi^)- Substitu- 
tion into (|1.8p yields the following Helmholtz eigenvalue problem: 

( A+(eA)^)*A = 0, r>l (3.4a) 

dr-^x = -iXPx, r = 1 (3.4b) 

-a*A = 37 + 4) (/3a, YS)YS - ^ (2 + ^s)Px , r = 1 (3.4c) 

outgoing r — s- 00. (3.4d) 

If A is such that p.4p has a non-trivial solution, then we call A a (deformation) 
scattering resonance energy or scattering frequency, and (4'a,/3a) a corresponding 
scattering resonance mode. 

Since outgoing solutions of the three-dimensional Helmholtz equation are linear 
combinations of solutions of the form (r)YJ'"(ri), \m\ < I, where h[^^ denotes the 
outgoing spherical Hankel function of order I, we seek solutions of the boundary value 
problem ([311) of the form: ^xir, fl) = A r,"(rj) h[^^ (eAr), ^xi^) = B Yf^fl), r > 1, 
G S*^, where A and B are constants to be determined. This Ansatz automatically 
solves the Helmholtz equation and satisfies the outgoing radiation condition. To 
impose the boundary conditions at r = 1 we substitute we substitute the expressions 
for ^A and /3a into p.4[) and obtain the following two linear homogeneous equations 
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for the constants A and B. This system has a non-trivial solution if and only if A 
solves the equation: 

ri e Xdh[^\eX) + X^h\^^ (eA) = 0, (3.5) 

where r; is given by (|l.lip . We shall be interested in the character of solutions to p.Sp 
for small, positive and fixed e. By (|D.11[) 

1- ^^S^ = 1- ..u, — = -(^ + 1) • (3.6) 

h^i^\z) r^+mpi{z) - pi+i{z)] 

Thus, we write p.Sp in the form: 

A^ + n G,ie\) = 0, G,(z) ^ f^J^iif). (3.7) 

h\ >{z) 

we call the solutions of (|3.7|) deformation resonances. Their properties are summarized 
in the following result, proved in Section |6l 

Theorem 3.2 (Deformation resonances). Fix e > and arbitrary. 

(i) There are 1+2 solutions of Equation (|3.7p denoted {A;j(e)}, fori = 0, 2, 3, . . . 
and j ~ 1, ...,/ + 2. These are the deformation resonance energies. 

(ii) The set of deformation resonance energies is a discrete subset of the lower 
half complex plane and is uniformly bounded away from the real axis. That is, for 
some h{e),K{e), 

5A/j(e)< 3Ai,j,(e) = SA^e) < 0, all / > 0, \j\<l + 2. (3.8) 

(iii) Aij (e) = 0{l) as I ^ oo. 

(iv) In the incompressible limit, e — > 0^, the imaginary parts of all resonances 
tend to —oo except for the family of Rayleigh resonances with real frequencies: 

A±(0) = ±y^i|^Ca+ 4(37-1), A±(0)=±y^4(? + 2)(? + l)(?-l), / > 2. (3.9) 

4. A theorem on resonance expansions and time-decay for the exte- 
rior Neumann problem for the wave equation. Our strategy for solving the 
initial-boundary value problem (jl.Sp . (|1.9p is to (1) construct the solution to the 
wave equation (|1.8ap with kinematic boundary condition (|1.8bp . which specifies Neu- 
mann boundary data on the unit sphere. This is the Neumann to Dirichlet map 
dtP ^ ^ = NtD [dtP]. Then, (2) substitute ^! = dtP NtD [dtp] into the dynamic 
boundary condition (jl.Scp and then study the closed nonlocal equation for /3 on the 
unit sphere. 

Denote by U the exterior of the unit sphere in M.^ : 

C/ = {x: |x|>l}. 

In particular, we consider the general initial-boundary value problem 

_ A) u 0, in [/ X (0, 00), (4.1a) 

u = 0, dtu = 0, on [/ X = 0}, (4.1b) 
d^u = dtf, on dU X (0, 00) (4.1c) 
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where the function dtf : dU x (0, oo) — > M is a prescribed boundary forcing func- 
tion, whose properties (anticipating those of the bubble shape perturbation, are 
prescribed below. 

We introduce a norm of functions, dtf ~ dtf{^,t), defined on the sphere, 5*^, 
which encodes smoothness in f2 G S*^ and decay in time, t. 
Definition 4.1. 

(i) For a function J : [0, oo) — > R, define for any k £ R. 

[JU = supe'^VWI- (4-2) 

t>o 

(ii) For a function 51 : Rt x S'^ — > R, with spherical harmonic expansion: 
g{t,n) - Y.i>o Y.\ra\<i9T{t)Yr{^), define, for g, k e R, the norm 

[5],.. - E E (1 + 0' ( isrwi + \dtgr{t)\ ) 

l>0 \m\<l *-° 
l>0 \m\<l 

The main result of this section is the following theorem, a generalization [ITIHS] . 

Theorem 4.2. Consider the Neumann initial-boundary value problem for the 
wave equation (j4.ip on the exterior of the unit sphere. 

(i) Time Decay Estimate: There exists a > 0, such that if [dtf]a,ri < oo, with 
77 > 0, then there exists a unique solution of the initial-boundary value problem equa- 
tions I^.Jp satisfying the bound 

' ' " " \0, \x\>l + ct 

Here, (Sw* < 0) denotes the scattering resonance of the exterior wave equation 
with Neumann boundary condition, which is closest to the real axis; see Theorem \3.1[ 

(ii) Resonance Expansion: Let a be as above. Suppose f : x R+ — > C has the 
spherical harmonic expansion in the sense that 

L 

hm [SJ - = 0, where f(^\n,t)^Yl E /rm™(f^)- (4-5) 

;=o \m\<l 

Denote by u^^\r,D,,t) the resonance expansion partial sum: 



u 



(i) 



-EE 



;=o \m\<l 



y I e-'i'-^^djr{s)ds — ("'-^:/^) 



k=l 



(4.6) 



Here h^^\z) is the outgoing spherical Hankel function of order I and 
{u!i} are the solutions of dhi^\iu/c) = 0; 

see Theorem \S.1[ The limit u{x,t) = liniL-i-oo w^^H^j ^) exists and converges to the 
unique of solution of the initial-boundary value problem (j4.ip . Furthermore, we have 
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the following error estimate for the resonance expansion. For < |a| + |^| < 2; 

^ fc[ai/-aj(^)]„,„^e---{''^i^"*|}(*-^), i<\x\<i + ct 

~ \0, |a;| > 1 + 

4.1. Proof of Theorem l4.2l . Wc first note that the decay estimate of Part 1 is a 
consequence of the resonance expansion of Part 2. Indeed, suppose Part 2 holds. Then, 
to prove the decay estimate (j4.4p . note that |u(a;, t)| < \u^^\x, t)\ + \u{x, t)—u^^^ {x, t)\. 
The first term is a finite sum, each of whose terms satisfies the decay estimate, while 
the second is controlled by (|4.7p . 

We now turn to the proof of the resonance expansion, Part 2 of Theorem l4.2l We 
first derive and, by Laplace transform methods, solve the equations for the spherical 
harmonic coefficients u™{r,t); see Equation (|4.24p . Then, we form the partial sum 
for it(^-'(r, ri, i) and prove convergence. In proving Theorem 14.21 we use the following 
two estimates: 

Proposition 4.3. Fix r > 1 and I >0. Then we have the bound 

- + (4.9) 

Proof. We use facts about hi{z), documented in Appendix ID. II Note that 
h[^\z) = z~^~^pi{z)e^^ , where pi{z) is a polynomial of degree /. Furthermore, we 
have the following recursion relating hi and dhi derivatives 

dh\'\z) = z-Hh\'\z)~h['^,iz); (4.10) 
see Appendix ID. 31 Therefore, 

pi{ujr/c) 



r'+i [lpi{uj/c) - pi+i{uj/c)] 



Now fix r > 1 and / > 0. The asymptotics (|4.9p follows using that pi{z) ^ a\z^ for 
large \z\ and that |a|/a|+}| = 1. □ 

Proposition 4.4. LetO< rj < =dof min;_fe |3a;;,fc|. Then, 

^ -1 r -1 

I — rj 



We now embark on the proof of Theorem 14.21 Substitution of the expansion 

00 

z.(r,r],i) = ^ ^ <(r,i)l^'»(17) (4.11) 
;=o |m|<i 

into equation (|4.1ap . expressing the Laplacian as + r^^A^ (radial and spherical 
parts) and using that — AsY;™ = l{l + 1)Y[''\ wc obtain the following equations for 
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ur{r,t) = 0. 



From the boundary conditions equation (j4.1c[) we have 

drur{r,t)^dtfr{t), atr^l 



(4.12) 



(4.13) 



In terms of the time Laplace transform of u, mI" 

initial-boundary value problem (j4.12l) . (j4.13p . with vanishing initial conditions for 
r > 1 becomes 



ipy + 

c 



'{r,p) = 



(4.14) 



druY^{r,p) ^ dtfl^ip), at r = 1, Mj™(r,p) as r oo. (4.15) 

Equation (|4.14p is solved in terms of the outgoing spherical hankel function /il^"*; see 
Appendix [D| 



ur{r,p)=Arip)h['\zpr/c), 



(4.16) 



where A"^{p) is to be determined. Imposing the boundary condition (|4.15p yields 



Arip) 



(ip/c) dh^j^^\ip/c) 



Inversion of the Laplace transform yields: 

1 r*"° „ajr{pW,'\ipr/c} 



Ai-ioo (ip/c) dh^p{ip/c) 



dp 



(4.17) 



where the contour is chosen in the right half plane, with JRp = /i > sufficiently large 
so that all poles lie to left. 

We find it convenient to work with a rotated (horizontal) contour and therefore 
make the change of variable 



Then, 



2n 



LO — ip . 

ddfi-zio)h['\u:r/c) ^ 

e auj 

(uj/c) dhl^^io/c) 



(4.18) 



(4.19) 



where Qi{uJ, r) is defined in (|4.8p . 

We first consider, r = |a;| fixed with r > 1 + ct and prove the finite propagation 
speed result: 



u™(r,t) 0, r>l + ct. 



(4.20) 
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Proof of Equation (|4.20p . By Theorem 13.11 the poles of Qi{ijJ, r) are in the lower 
half u— plane. Therefore, for any p > 



+ I ] e-*'"(*-('''i)/^) Q/(w,r)dw = 0, a* > 0, (4.21) 



1 

where Cp ~ { uj : \uj — ip\ — p, 3w > yU }. It follows from ()4.2ip that for r > 1 + c< 

u^{r,t)= lim — / e'"«'-i)/= - *'Q,(tj,r)dtj, 

P-oo27r7_,+,, (4.22) 

= - lim — / e*""'-!)/^ - *)Q,(w,r)dw 
p^oo 27r 

Note that the integrand is analytic in the upper half plane, continuous up to the real 
line and, by Proposition 14. 3[ maxgg[o,7r] IQKp^*^, ?')| — ?> as p -H' oo. By Jordan's 
Lemma the limit in (|4.22p vanishes, proving (|4.20p . □ 



For r < 1 + ct, we decompose the Laplace transform of the initial data, dtf. 



rt-{r-l)/c 

9J-(-jc.)= / e^'-^dsfr{s)ds+ I e'^'dsfr{s)ds (4.23) 

t~(r-l)/c 



and the corresponding solution 

J-oo+ip {u/c)dh\ \u/c) 



1 



(4.24) 



2^^J-oo+^^^ ■ (CU/C) dhY\L0/c) 
V ' 

In the following two propositions, we evaluate ?7/|^| and J//^^. 

Proposition 4.5. Let ui^k = ^i,k (c~^) , fc = 1, 2, ...,/ + 1 denote the solutions 

of S/ip' (oj/c) = 0; see Theorem \ 3.1\ with e = c^^ Q. for r < 1 + ct, Ui^m can be 
expressed as a sum over residues: 



fe=i -'o "^-^f (w/c) a/ij (w/c) 

. - r"-"%,-..;..-.,a,;™„„, '■!"'-".'[/" . (4,25) 



''The notation of Theorem 13.11 suggests we label the solutions of dh^^\u)/c) = as uii ^ (c~^). 
For ease of presentation, in this self-contained section, we mildly abuse notation here and refer to 
these solutions as simply ui^k- 
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Proposition 4.6. Ui'^^{r,t) = 



(1) 



Proof of Proposition \4-5\ We evaluate the contour integral representation of ^ 
by the method of residues. Consider the clockwise-traversed rectangular contour 
pM,M,7 = Ti + r2 + + where 



Ti = {x + iii: -M <x<M}, = {M - iy : -fi < y < 7} 
Ts = {-X - ij : -M <x< M}, T4 = {-M + iy : < y < fi}. 

By Theorem EH dh[^\u}/c) has l + l zeros in the lower half plane and we can therefore 
choose M > and 7 > such that r^^'^-'' encircles these zeros. 
Claim 1: Jp^ and Jp^ — > as A/ — > 00 

This claim is a direct application of the asymptotics (|4.9p on the Hankel function ratio 
in the integrand. Namely, 



t-(r-l)/c 



r-i h\^\ujr/c) e 
" as 



< 



t-(r~l)/c 



e'^(*-^)|9,/,"| dse- 



(uj/c) dh^\uj/c) 
>,.,iizii 2 



■ duj 



\ujr/c\ 



du) 



< 



C 

rM 



i-(r-l)/c 



\dsfr\ ds ^ aa M ^ 00. 



We therefore have 

l+l /.t-(r-l)/ 



l+l rt-(r-l)/c , (1) / I \ 

Jo (yijc) dh) ' (uj 



k=l 

1 



27r /rJ^-f 



t-(r-l)/c 



e-(«-^)g./r(.)d. , (4.26) 

(w/c) araj ' (w/c) 



where M > sufficiently large. The size of M depends on I, since the contour must 
enclose alH + 1 poles of the integrand. This completes the proof of Claim 1. 

Our goal is to deform the contour T3 = T^^''' downward, by sending 7 ^ 00, and 



to establish the following claim, from which Proposition 14.51 follows immediately. 
Claim 2: The doj integral in (|4.26p satisfies the estimate 



sup 

M>0 



duj 



0, as 7 — > cxj. 



(4.27) 



To prove Claim 2, we integrate by parts in s, exploiting the oscillatory character 
of the integrand, in order to obtain absolutely convergent duj integrals. 



• • du) 



1 

27ri 



[hir.t) + hir,t) + h{r,t)] 



(4.28) 
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M-i7 pt-{r-l)/c 



where 
h{r,t) 

-M-i-y Jo 

I,ir,t)=-djrit-ir-l)/c) 



'''-'dlfr{s)ds '\^':^'^ ^^ <^ (4.29a) 
(oj^'/c) ah] ' {bj/c) 



duj (4.29b) 



■ duj 



(4.29c) 



We conclude the proof of Claim 2, and therewith Proposition l4.5i by now showing 
forj = l,2,3, sup|/j(r,t)| = 0(7"^) ^ 0, as 7 00. (4.30) 

j\/>0 

Estimation of Ii{i\t): Choose 7 > 2r]. Using the bound (j4.8p . we have 



\h{r,t)\< 



Al~if pt-(r-l)/c 



< 



M-i^ Jo 
t-(r-l)/c 



^^/i[%;r/c)e 



-— 



uj/c)dh[^'> (uj/c) w 



2c 







dw 



8c 



7r 



Estimation of l2{r,t): 
\h{r,t)\<\djrit-ir-l)/c)\ 



— duj 



M-i7 I'^r 



{u:/c)dhf^ {uj/c) 

ryj- 



Estimation of {r, t) : Lastly, 

f.A/-i7 



\l3{r,t)\<\dtfr{0)\ 



< - mr 

r 



g-iwtgiw((r-l)/c) 



lit-{r-l)/c) 



M—i-f 



[i0/c)dh\'\L0/c) 



diu 



dw 



< 



2c 

7 r 



-7(t-(r-l)/c) F'^ 



Proof of Proposition \4-()\ We must prove that Ui^l^{r,t) = 0. Note the following 
bound, for 5a; > and r < 1 + ci: 



e'^'dsfr{s)ds 

t-{r-l)/c 



< C|c^rie-(^'"+")[*"(''"^'/"l[9t/]a,^. (4.31) 



This follows from direct substitution and integration of the bound 

|e*--a,jr(s)| < e-^-^"^ e-^'^ [dtfU^. 
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We now use (|OT|) to establish that J7/_^(r,t) = 0. Indeed, 



1 

2^ 



M+if 



h'l^\ujr/c) 



< 



2Cc 
27rr 

< C 



-M+iu 
M+iv 



(u;/c)5/iW(^/c) 

-{v+Ti)[t-{r^l)/c\ 



duj 



-M+iu 



M4 



-r]{t-(r-l)/c) ^ fj 



duj 

°° dz 



dz 



(2'] 

as V ^ oo, and thus, Ui^{r,t) = 0. This completes the proof of Claim 2 and thus 
Proposition 14.51 □ 

Proof of the Resonance Expansion ()4.5p . The resonance partial sum, defined as 



L 

EE 

Z=l |m|<i 



.£t^o -^f (../c) 



which gives Equation (j4.6p . 

To study the convergence of m'-^' as L tends to infinity, consider the difference 
y(^) _ u'^^\ where N and L arc large with N > L and prove: 

Proposition 4.7. Let |3a;*| = mmk.i \'^i^i,k\ and pick i-j £ (0, |3cli*|). Assume 
[dtfU^ < oo. Then, for \a\ + \b\ < 2, 



(4.32) 



Therefore, for any compact subset Itxflxd x R^, {m'^^} is a Cauchy sequence in 
C^{It X fix) and converges uniformly on If xft^ to a limit u(x,t), which is a classical 
solution of the initial-boundary value problem (|4.1I) . 

Proof. We begin by estimating the difference u'^'^\x,t) — u'^^\x,t): 



u'^^\x,t)~u^^\x,t) 



N+1 

< E E 

l=L+l \m\<l 
N+l l+l 

< E E E 

l=L+l \m\<l 



E 

k=l • 



k=l 



e-'"' ^'-'^dsfris) ds 

< \Yrm\\dfrL.,T. I 



fc=l 



(4.33) 



To complete the proof of Proposition 14. 71 we require Proposition and the following 
two propositions: 
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Proposition 4.8. There is a constant C > 0, independent of r, such that for 
< a < 2, and any r > 1, 



< C 



(1 + ly* ■ l" 



(4.34) 



Proposition 4.9. |9^9^Y(™(6i»| = 0(/i/2+n)^ for\a\ + \h\=n, n = 0,1,2. 

Proposition 14.81 is proved using the uniform asymptotic expansions for Hankel 
of large argument and large order and the asymptotic locations of the resonances 
{uji^k}- The arguments are of the type used in [5T] and very close to our proof of 
Proposition 15. 2[ presented in Appendix [F] Proposition 14.91 is a classical result of 
Calderon and Zygmund [5] . We now complete the proof using Propositions 14. 4[ 14.81 
and 14.91 First, we bound the term of ()4.33p as follows: 



\yrm\\dfr\^,, E 



fc=i 



h^p (uji^kr/c) 



(1 + 0^* 



Summing this estimate over / and m yields: 



{x,t) 



< 



C 



[dtf 



(N) 



dtf 



Similarly, we have for any a, b with |a| + |6| < 2, we obtain (|4.32p . This completes the 
proof of Proposition 14.71 Theorem 14.21 now follows by passing to the limit — )• oo, 
which proves the resonance expansion bound on u — u^^\ (|4.7p . □ 

5. Linearized bubble dynamics: solution of the initial-boundary value 
problem, decay estimates & resonance expansion. We next state our main re- 
sults on decay estimates and resonance expansions for solutions of the initial-boundary 
value problem (|1.8|) . p. 91) . A key tool is the general result, Theorem l4.2[ for the Neu- 
mann to Dirichlct map. 

Theorem 5.1. Fix e > and arbitrary. Suppose \\/3{t = 0)|| = J2i>o J2\m\<ii^ + 
1)^+6 1/3™ (0)1 < oo. There exists a unique solution 4'(r, il,t), /3(f2,t), defined for 
r > 1, ft G , which solves the initial value problem (|1.8p . 

(i) Decay Estimates: The solution satisfies the following decay estimates: 

1/3(1], i)| < C ||/3(t = 0)11 e-i^^*(^)i*, nes^, (5.1) 
'C ^e-™{^^*('^)>l'^'^-(<^)l>(*-^(l^l-i) ||/3(i = 0)||, l<\x\<l + e'H, 



\^ix,t)\ < 



0, > 1 -t-e^H. 

(5.2) 



Here, (e) < and (e) < are, respectively, imaginary parts of the deformation 
and rigid scattering resonances, which are closest to the real axis; see Theorems \S.1\ 
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(ii) Resonance Expansion of /3: Define the resonance partial sum for the 
bubble surface perturbation: 



L 



l=a \rn\<l 
L 1+2 

= /3r(0)E A,,e-^^'-* Res [ X'+riG^eX) Yrm, 



X=Xi 

1=0 \m\<l 3 = 1 



(5.3) 



where {Xij}, the solutions of +riGi (eA) = 0, are the deformation scattering res- 
onances; see p.Sp and Theorem \3.'A The limit [i[VL^t) = Mitil^oo 1^^^^ {^,t) exists 
and converges uniformly in (7^(5*^). 

(iii) Resonance Expansion of ^P: Given (3, defined in (ii), define 



1=0 \m\<l k=l 



[e^uJi^k)d'^hf'\eu!iM) 
(5.4) 

where {uJi k\ denotes the set of rigid Neumann scattering resonances; see Theorem lS. 11 
The limit ^'(r, f2,t) = lim^^cc ^''■^■'(j', £7, t) exists and converges in G^ on any com- 
pact subset _ftr C {|x| > 1}. 

Proof of Theorem \5.1\ We apply tlie results of Section 2] with u(r, f7,t) = 
*(r,rj,i), f{Vt,t) = l3{Vt,t) and c = e-\ with dr'i!{l,Vt,t) = dtl3{n,t) as Neumann 
data, coming from the linearized kinematic boundary condition and obtain: 

oo 

^ir,n,t)= J2 E *™(^'^)^"(^)' (5.5) 

1=0, 1^1 \m\<l 

-, pti-\-ioo h^^Uifnr) 

^i,^{n,r,t) ^ — eP'dtPriP) ' m / = 0,2,3,..., \m\ < I 

^■^iJp.-ioo [iep) dh\ ' {ie) 

(5.6) 



Setting uj ^ip in (|5.6p and applying Proposition 14.61 we obtain 

^'J"(17,r,i) = — / e"*"* / e*-^" 9,/3["(s) ds ^ , ^ duj . (5.7) 

27rJ„^+j^ 7q Gi{euj) 

Our goal is to obtain a closed equation for /3. This we obtain by substitution 
of (|5.5I) , (|5.7|) into the linearized dynamic boundary condition (|1.8cP . Evaluating (|5.7p 
at r = 1 and expressing it as a time convolution, we get, for each / = 0, 2, 3, . . . and 
|m| < /, 

t / 1 



dt^rih^,t) ^ dt ^C-'-^^j{t-s)d,pr{s)ds, (5.8) 

where £~"'^ denotes the inverse Laplace transform. Now taking the Laplace trans- 
form of equation ()5.8p we obtain for the Laplace transform of the linearized dynamic 
boundary condition (|1.8cP : 

^1—^ [pWHp) - P/3r(0)] = [37 + ^) Sio + ^(/ + 2)(/ - 1)] ^"(p), (5.9) 
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for / = 0,2,3,..., and |m| < /. In obtaining (|5.9|) . we have used well-known rela- 
tions for the Laplace transform of a time-derivative and a time-convolution, and that 
-AsY," = 1(1 + 1)Y,". Next, solving ([SI]) for /3p(p) yields 

^(p) = /3r(0)p {p' - riGi{tep)y\ (5.10) 

where n = Wci^ + — 1) + 87 + <5;o. Inverting the Laplace transform and 
changing variables to A = ip, we have 

Pnt) = /3r(0) - / e-'^* dA, /z > 0. (5.11) 

27rJ_^+j^ A^+nG/(eA) 

Thus, 

mO^EE^^y + .>0. (5.12) 

We will show (1) each term in the series (|5.12p can be expressed as a finite sum of 
residues at poles in the lower half A-plane and (2) under suitable regularity hypotheses 
on /3(ri, 0), the series converges uniformly in C^{Sq x M^). 

Evaluation of /3™ (t) , in terms of residues: The poles of the integrand of (j5.1ip 
are precisely the deformation resonances, described in Theorem 13.21 In particular, 
there are I + 2 poles, zeros of A^ -I- riGi{eX). 

Choose r*^'^'''' = Fi -I- r2 + Fa + r4 a rectangular contour of the type used in the 
proof of Proposition 14. 5| with M, 7, fi so that F*^'^''' encircles these I + 2 poles in the 
complex A-plane. By the residue theorem 

VzA,,e-^^'-* Res [X^+nG^ (eA)]"' = --^ / "T^^* d\. (5.13) 



Note that 



A"(0=/3r(0) lim / ■■■ dX, (5.14) 

where • • • denotes the integrand in (|5.13p . Thus, we would have 

1+2 

PTit) = V A,,,e-'^'-* Res [ A^ + nGi (eA) ] , (5.15) 

0=1 

as well as the expansion (|5.3p , if we can prove /p^ ^ - ■ ■ dX — > as M — )■ 00 , and 
Jp • • • dA — >■ as 7 — > 00 . 

To prove /p^ ^ • • • c?A — > as Af — > 00 we begin by noting (Appendix ID.ip the 
identity 

Gi{eX)^-{l + l)+ieX+'-^^^. (5.16) 

Thus Gi(eA) to be asymptotically linear in eA. Thus AI can be chosen large enough 
so that 

^ < TTTJ, for ^X > M. (5.17) 



|A2+r,G/(eA)| " |A|2^ 
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Therefore, we have 



A2 + nGz(eA) 



dX 



< I |e-'^*|^|^dA<G(Ai+7)e^*^^0 as M^oo. 



A similar estimate holds for the integral along 

We next bound the integral: limA/_^oo /p^ ■ • ■ dX = /cj_)^_ ^ ■ ■ ■ dA as 7 — >■ oo. For 

7 = 3A sufficiently large, 



A + 



riGijeX) 
X 



1 



-{l + l) + ieX + 



e\p\{eX) 



(5.18) 



Al + enC'(A-i) A 



[l + enO(A-i)] 



Therefore, dx \ X 



riGijcX) 
A 

tegration by parts we have 



-A 2 + C'(A ■^). Using these asymptotics and in- 



iXt 



iX 



M-t-y X'^+riGi{eX) 



< 



< 



1 

t 

nie 



-i\t 



t 



dX 



--+0{X~')]dX 



i(l + .n0(A-))irA 



■dX + 0{M-^) < 



me 



'It 



-, i>0. 



(5.19) 

Recall that 

r j-M-ia -1, 

/ ■■■ dX= \im e-*^*-^ TTTTT'^^- 

JaA=-7 ^'^°°J-M-to^ X^+riGi[eX) 

Thus, passing to the limit in (j5.19p . first M — > 00 and then 7 — > 00, we obtain that the 
contribution from the contour can be made arbitrarily small. This proves (|5.15p . 
the expression for (3p{t) as a sum over residues. Furthermore, multiplication of ex- 
pression (|5.15p by Yi"^{il) and summing over I from up to L yields the expression for 
/?(^)(f],t), the resonance expansion partial sum for the bubble surface perturbation. 

We now discuss the convergence of P'^^^ (fi, t) as L — >■ 00. Let N > L and consider 
the difference: 

AT 1+2 . 

/?W(i7,t)-/3(^)(i7,t) = E E prmY^e-^'-^' Res rrm. 

l=L\m\<l j=l + ; il J 

(5.20) 
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We require the following bound, proved in Appendix IF] 
Proposition 5.2. Fori sufficiently large, 



Res — = O (l-^") . (5.21) 



To bound the difference in (|5.20p we use Proposition 15.21 and Theorem 13.21 char- 
acterizing and estimating the deformation resonances, Ajj- and the pointwise bounds 
on Y[^{U) of Proposition 14.91 Let denote any differential operator of order 
acting tangentially to 5*^, i.e. 5^ = 9g^d^^, \a\ = ai + a2- Then, 

N 1+2 

l=L\m\<l j = l 

N 

l=L \m\<l 

(5.22) 

Thus, if E^lLE|m|<i(l + 0^+^ l/3r(0)l < we have that /3'-^Hn,t) and 
derivatives in and t up to order two, converge uniformly in Sq x M.f . Moreover, the 
liuiL^oa Pi^,t) is in C^^S^^ x M^) and satisfies the decay estimate 

N 

I P{n,t) I < E E (1 + 0'+^ Wm ■ . (5.23) 

l=L \m\<l 

Now apply our result on the Neumann to Dirichlet map, Theorem l4.2l . to conclude 
for |x| > 1, t> Q that 

< ^^-.ain{m,ie)\}t ^ ^ (1 + /)2+i |;3r(0)| . ^^'^^^ 
;=L |m|<i 

The second inequality in (|5.24p . follows since |3A*(e)| < |5a;*(e)| = 0{e~^) for e 
sufficiently small. Moreover, (\l>(r, f2, t), /3(51, t)) is a classical solution of the initial 
value problem. 

6. Scattering Resonance Energies — Statements of Detailed Results. 

As discussed in Section [Sj two families of scattering resonance energies arise in the 
analysis of the linearized bubble dynamics near the spherical equilibrium state: 

1. Rigid resonance energies, {wi^fe}, solutions of 9/ip^(ea;) = 0. 

2. Deformation resonance energies, {Xi,k}, solutions of 

z dh'-^Uz) 

X^+riGiieX)=0, Gi{z) = ^if^, n is given by dm]). (6.1) 

h\ (z) 

Both classes of resonances depend on e. The incompressible limit, e — > 0, is a 
singular limit. Indeed our results show that for each ? = 0, 2, 3, . . . , as e tends to zero, 
the imaginary parts of all Z + 1 rigid resonances and of I of the I + 2 deformation 



5^(/3W(f),i)"/3^'^)(f^,t)) 
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resonances lie near circular arcs in the lower half plane and have imaginary parts 
which tend to minus infinity. The remaining two deformation resonances are just 
below and converge to the real axis as e — ?► 0; see Figure [Ol of the Introduction. The 
imaginary parts are 0{e~^ exp(Ke~^)) as e tends to zero. 

6.1. Deformation resonances; {A;j(e)}, such that + r; G/(eA) = 0. 
Theorem 6.1. Fix e > and arbitrary. 
(i) There are I + 2 solutions of 



(6.2) 



denoted {Xi j{e)}, for Z = 0, 2, 3, 4, . . . and j ^ 1, ... ,1 + 2. These are the deformation 
resonance energies. 

(ii) {A/j (e)} are symmetric about the imaginary axis and satisfy 3A;.j < 0. 

Theorem 6.2 (Asymptotics of Rayleigh deformation resonances, A;*(e) for small 
Mach number, e). 

(i) For I ~ 0, the two deformation resonances are located in the lower half 
plane, a distance 0{e) from the real axis, are given by: 



Aif(e)=±Wro 



ro 2 •'"0 
4 2 ' 



ro = |^Ca+2(37-l);^ . (6.3) 



(ii) Fix > 2. There exists e*{L*) such that for e < e* and 2 < I < L* , two 
of the 1 + 2 deformation resonances (see Theorem \6.1\) . are located very slightly below 
the real axis and given by: \f — 5RA^(e) + i'^Xfie), where 



mf{e) = ±J—il + 2){l-l){l + l) 



We 



a+2)(/-i) 

2{2l - 1) 



9A±(6) = -- 



-[(/ + 2)(/-l)]'+i(? + l)' 



2H\ 



(20! 



'We. 



VWe 

21+2 



(6.4) 



(iii) In the incompressible limit, e — > 0+, the deformation resonance equa- 
tion (j6.2l) . which has I + 2 roots in the lower half plane for e > 0, reduces to the 
quadratic equation: 

X^+n GKO) = A2 - 37 + ^) 5i^-^^{l + 2){l + l){l-\) = 0, Z = 0,2,3,... 
and has two real roots: 



Ca 2 
2" ■ We 



1 

We" 



Af (0) = ± J37 ( ^ + ^ ) ^io + ^(/ + 2)(/ + 1)(? - 1), Z = 0, 2, 3, . . . (6.5) 



which are real frequencies. The corresponding solutions of the linearized time- 
dependent perturbation equations: 



are undamped and oscillatory solutions. 



± g-»A±(0)t y^m(^)^ /=0,2,... 
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Remark 6.1. In Section^ we present an asymptotic analysis and numerical 
computations offering strong evidence for the scattering resonance energy nearest to 
the real axis being a distance of order 

^(e-i Wee-^e"'"^""") , k > 0, occurring at order /^(e) = ©(Wc e'^) . (6.6) 



Finally, the following theorem establishes that for any fixed e and sufficiently large 
I, the deformation resonances (excluding one) are distributed near a circular arc. In 
particular, this implies that for any fixed e scattering resonances are in the lower half 
plane and uniformly bounded away from the real axis. 

Theorem 6.3. Fix e > and arbitrary. The I + 2 deformation resonances are 
approximated as follows: 

(i) There are I + 1 resonances near zeros of dH^\e\); see Equation {E.IT^ . 
Precisely, there exists Ki > such that for I > Ki (e~^We)^^^^~^'^ , < q < 1/3, 
those in the fourth quadrant can be expressed as: 

1 + 1/2 



l + 2-i/3e-2"/3(; + i/2)-2/3|^,| 



+ ^2-i/3e-4-/3(; + i/2)-4/3 + o(^5-i/3(; + i/2)-2/3-9) 



•7) 



for s = 1,2, ... , [{I + 1)/2J + 1. Here, rj'^ is the sth zero of dzAi{z); > 77^ > 773 > 
• • • J22^ . Note the error is uniform in e. These resonances come in pairs, symmetric 
about the imaginary axis, {Xls, —^Ls}- For even I, Xi^(i+2)/2 *s « simple resonance on 
the negative imaginary axis. 

ill) The remaining deformation resonance is located on the negative imaginary 
axis and can be expressed, for some K2 sufficiently large and I > K2 e~^We, as 



^1,1+2 



(6) = ^ 



.1 + 1/2 



e^{l + l/2) 1 We 
m ^ 2? 



^+0([.^Wc-'(, + l/2,]-) 



(6-8) 

(iii) The set of deformation resonance energies is uniformly bounded away from 
That is, for some Z*(e),j*(e), 



3Ai,j(e) < 3A;^j,(e) ee SA^e) < 0, all I > 0, \j\ <l + 2 



(6.9) 



(iv) Moreover, there exists C > such that as I 00, A;j(e) = e ^0{l), 
C e-i?i/3 < |3A,,,(e)|. 

We conclude this section with a proof of Theorem 16.11 Theorem 16.21 is proved 
in Section |8l Theorem 16.31 is proved in Section [9l To prove Theorem 16.11 we use the 
following rational function representation of Gi{z) = dh^^\z) / h^^\z) , which follows 
from (|D.5I) : 

Proposition 6.4. Gi{z) = -(/ + l) + iz + ^^f^. 
Proof of Theorem \6.1\ 
(i) We first observe that there I + 2 scattering resonance energies. To see this 
we observe, by Proposition l6.4[ that A'^ + riGi{eX) = is equivalent to 



ri 



-{I + 1) +ieX + 



Pii^X) 



= 0, 



(6.10) 
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which clearly has I + 2 roots since pi{z) is a polynomial of degree I. For e = 0, the 
equation reduces to A'^ = r;(Z + 1), yielding two real frequencies, corresponding to the 
non-decaying Rayleigh modes of the incompressible limit problem |13| . 

(ii) We now prove that if A is a solution of (|6.10p . then so is —A. By ()D.4[I 

L eXdHl^\e\) 



X' + n Gi{tX) = ^ X' + n 



Hi'\eX) 



0, 1 + 1/2. (6.11) 



Using (|D.10|) we have A^ + r; 



0. 



Taking the complex conjugate and using the analytic continuation formulae (jD.7[) 
and lESl) we obtain 



= A + r; 



61/ — eA 



Hi'\eX) 



n 



1^ ^ ( t A+m-^)) 

- + eiy + e(-X) — 7-T ^— 

2 Hl'\e{-X)) 



(6.12) 



Thus, —A satisfies (|6.1ip and is therefore a scattering resonance. 

We now prove that all resonances lie in the open lower half plane, 3A < 0. We 
begin, following [IT], by showing that there are no resonances along the real axis. Let 
ly = I + 1/2. Using the definition of H^^^ in terms of Bessel funcitons, equation (|6.1ip 
for the resonances can be written 

= (^^ - [J, (eA) + lY, (eX)] + (eA) [dJ, (eA) + idY, (eA)] . (6.13) 

Now assume X = x/e, where x is real and nonzero, is a resonance and we show that 
this leads to a contradiction. Indeed, 



1 1' 

^ t;] [J-^i^) + iY^{x)] + X [Jl{x) + iY^{x)] 



£2 n 



£2 n 





Y 1 x^ 




+ i 


\£2 n 


-I) 



Since x G M both quantities in braces equal zero. Therefore, we have 

1 x^ 1\ xJlix) xY/ix) 



2 

£ ri 



Mx) 



Y,(x) 



(6.14) 
(6.15) 



implying —x [ Jl(x)Y^(x) — Y/{x)Ji,{x) ] ~ 0. The expression in square bracket is 
the Wronskian of two linearly independent solutions of Bessel's equation and satisfies 
the following identity [1]: 



J'M'y^x) - Y/{x)J,{x) = — ^ 0. 



(6.16) 



Hence, for a; 7^ 0, we have a contradiction and thus 5A ^ 0. 

Next, we claim that there are no resonances in the upper half plane, and therefore 
5A < 0. For if A is a scattering resonance with 3A > 0, this would give rise to a 
solution of the initial-boundary value problem for the linearized bubble perturbation 
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system (Ull) of the form: *(r,17) = e~'^*h\^\eXir), f3{n) = e-'^'y™(17). Since 
3A > 0, this solution would decay to zero exponentially as r = — > oo ( is 
e'^ times a rational function of z), be smooth and be exponentially growing in t. This 
contradicts conservation of energy, Proposition 12.11 Thus, 3A < and by the above 
discussion 3A < 0. □ 

7. Longest-lived mode; formal asymptotics and precision numerics. 

From Theorems 16.11 to 16.31 we know that for any finite (not necessarily small) e, 
there is a strip: 

5A±(e) < 3A < 0, 

in which there are no scattering resonances. Moreover, the incompressible limit spec- 
tral problem, e = 0, has eigenvalues only at the real Rayleigh frequency pairs, given 
in Equations (|6.3p and (|6.5p . and none in the lower half plane. 

Question: How thin is the strip? That is, can we estimate 3Aj(e), for e small? 

Figure [TTl] displays the distribution of numerically computed Rayleigh resonances 
for a fixed e > and several values of I. The figure shows the resonance in the 
lower half plane (a pair, Aj(e) by A i-^- —A symmetry), which is closest to the real 
axis. These resonances correspond to the slowest time-decaying term in the resonance 
expansion and the resulting decay rate in Theorem 14.21 

Let U{e) denote the mode index of the resonances nearest the real axis and 
\f{e) = {e) ^^'^ resonance itself. In Equation (|6.4p we display the leading order 
behavior of the QA;*(e), for / not too large, and e small: 

9Af(e) iM{l,e) (7.1) 

where 

Mil,e)^^Jil + 2)il-l)t'il + iy 

We attempt to compute the h{e) and !3Aj(e) ~ —A4 (Zj,(e),e), by minimization of 
A4{l,e) over I, treated as a continuous variable, while keeping eWe~^^^ fixed. The 
calculation presented below in Remark 17.21 suggests the following conjecture: 
Longest lived mode index: There exists A, > such that for e < e* 

h{e) ^ Ae-^We. (7.3) 

Real part and imaginary part ("reciprocal lifetime") of longest lived resonance mode: 
There exists ri,ii,i2 > such that Xf{e) = ?R.Xf{e) + i3Aj(e), where 

5i{^?(e)}~±^(^"'We), 5{A±(e)} « (e-2we) exp(-^2e-'We) (7.4) 

Remark 7.1. |3Aj(e)|, the distance from Aj(e) to the real axis, tends to zero 
very rapidly (beyond all orders in e) as e ^ 0. In contrast, the radial resonance, 
Ao(e), tends to the real axis only linearly in e; see Theorem \6.2[ 



2H\ 



/We 



(7.2) 
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Numerical evidence supporting (|7.3p and 17.41 is obtained by numerically com- 
puting the resonances for fixed e and various I, and determining the optimal mode 
index. Z+(e), and resonance 

)^numericai(^^y The Mathcmatica function FindRoot with 
options AccuracyGoal -> 20, PrecisionGoal -> 20, WorkingPrecision -> 2000, 
and Maxlterations -> 1500 were used to compute the resonances. Figure mi shows 
the numerically computed optimal mode along with a numerical fit of the data. Note 
the expressions in Equations ()7.3|) and ()7.4p are the leading terms of the numerical fit 
functions. 




Fig. 7.1. Numerical computations of (a) angular momentum index of slowest decaying 

mode (left column), (b) real part, 5RA*(e) (center) and (c) imaginary part, QA*(e) (right) of the 
corresponding Rayleigh - deformation resonances, ^i^(e) = ^*{e). Numerically computed values are 
indicated by dots. A numerical fit (solid curve) to the functional form displayed in the inset of each 
plot is also given. 



Remark 7.2. Formally minimizing A^(Z,e) over I > 2 for fixed eWe"^, will 
approximately determine the longest lived mode and its lifetime. The optimization 
can be performed by treating I as a continuous index and applying Stirling 's formula, 
71 ! « \/2nn (n/e)", n ^ 1. This leads to a formal approximation of /*(e), the 
angular momentum index of the most slowly decaying mode, and A^^^j = = 
3?Aj(e) -|-j9A^(e), < 0, the corresponding scattering resonance energies. In 

particular, we obtain: h{e) - ^ e^^ We, 3fiA±(e) - ±i {^f^^ (e"^ We), 3A±(e) ~ 



8. Rayleigh deformation resonance analysis; Proof of Theorem 16.21 In 

this section we prove the assertions of Theorem 16.21 concerning the Rayleigh reso- 
nances, A;^(e), Z = 0, 2, 3, . . . appearing very close to the real axis for e small. 

The main ingredients of the proof of Theorem l6.2l are the following several proposi- 
tions, which establish the detailed structure of the Taylor expansion of X^ + riGi (eA) = 
for A bounded and e sufficiently small. This will enable the explicit determination 
of the Taylor expansion of the Rayleigh resonances, A^(e), and in particular their 
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imaginary parts. 

Proposition 8.1. For e small, the equation + riGi{e\) = can be written as 



^ , ^ ^ ' 



= 0, (8.1) 



where A{X) is even and real for real X and B{X) is odd and purely imaginary for real 
X. 

Proof. This proposition follows by Taylor expansion of Gi{z) about z = and the 
following: 

Claim 8.2. For all j = 0, 1, 2, . . . 

5|gP^(0)} = 0, 5R{gP+'^(0)} = 0. (8.2) 

This claim follows directly from the observation that 

Gi{z) is real on the imaginary axis. (8-3) 

To prove (|8.3p . we consider Gi{iy), y G K. By Proposition [531 it suffices to show that 
the ratio p'i{iy) / Pi{iy) is purely imaginary. By (|D.6[) we have 

Piiw) = XI °^n{iyT = XI 2'-" a -nMn! ^"' ^^^^^ purely imaginary 
The ratio is therefore purely imaginary, which implies that Gi(iy) is real for y real. 



We now prove Claim [8?2] bv induction. Consider the rational representation of Gi 
in Proposition 16.41 As noted in p.6p . G;(0) = —{I + 1). Moreover, for y S M 

dGm = li- ^iM^:^ . nm ^iM^^ , 

y^O iy J/-S-0 y 

since Gi{iy) € E for y € E. 

Next, we make the induction hypothesis: 

» (^^(0) j - 0' (^^27^(0)] = (8.5) 

and seek to prove the assertion (|8.5p with j replaced by j + 1. 
From the Taylor series of Gi{z) we have: 



^^^^'-(0) ^ lim 



^^(2fc-2)! dz2'=-2 ' ' ^(2/c-l)! (iz2fc- 

(8.6) 

Thus, 3 d^^Gi{Q) = is real by induction hypothesis 
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Similarly, we have 



T^(0) = lim 



(2j + 1)! 



^ (2/t)! dz2'= ' ' ^(2fc-l)! dz2fe- 



(8.7) 



which is purely imaginary by induction hypothesis and since 5 d^-'Gi{0) = 0. This 
completes the proof of Claim [5^ and thus Proposition 18. II □ 

Proposition 8.3. There is a positive integer 2 < < oo, for which Equa- 
tion 118.1]) is of the form: 



G 



(2g.+i; 



(0) + O {e^^'+^) 



0, 



■A(A) 



e(A) 



where A{X) is even and real-valued for real A and B{X) is odd and imaginary-valued 
for real X. 

Proof of Provosition [87^ We begin by observing that G;(z) is not an even fimction 
of z. This is a direct consequence of the asymptotic behavior of G/(x) as x ^ ±oo: 



Gi{x) ^ -{1 + 1) + ix + 1 = -1 + ix, a;>l 
Gi{-x) ^ -{I + 1) - ix - I = -{21 + 1) - ix, X > 1. 



(8.9) 



Therefore, there is a first non-zero odd derivative, Taylor coefficient. We call the order 
of this coefficient, q^,. The next result asserts that q^, = I. □ 
Proposition 8.4. 



Gp^+'^(0) 



z (2/ + l)!(^)', k = l, 

0, 0<k<l-l. 



(8.10) 



Proof of Proposition \8.4\ Since we are interested in the odd derivatives of Gi at 
the origin we calculate: 

Giiz) - Gii-z) ^-il + l) + iz + + + + zMtil 

Pi{z) Pi{-z) 



2iz + z 



Pi{z) Pi{-z) 
2ipi{z)pi{-z) +p[{z)pi{-z) +p'i{-z)pi{z) 



Pi{z)pi{-z) 



111) 



By jESl), a\ = i-'-i and therefore 2ipi[z)pi{- z) = 2ia\z a\{-zy + ©(z^'-i) 
—2iz'^^ + 0(2^'^^). This suggests separating out this term and writing 



Giiz)-Gii-z)^ z 



2iz2' + Ri{z) 



Pi{z)pi{-z) 



.12) 



where 



Ri{z) = 2ipi{z)pi{~z) + p[{z)pi{-z)+p[{~z)pi{z) - {-2iz^') 



.13) 
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The expression (|8.12p simplifies due to a subtle cancellation which follows from 
Lemma 8.5. 



Ri{z)=0 



(8.14) 



and therefore 



Gi{z)-Gii-z) 



-2iz 



21 + 1 



_pi{z)pi{-z) 

We prove (|8.14p below, but first note its consequences. For z small we have 



Gi{z)-Gi{~z)^ 



where we have used that 



Since 



-2iz^'+' 



r oi 



A) +oiz) 



= 2i 



2H\ 



{21)1 

(20! 
2H\ ■ 



^21 + 1 



(8.15) 



Gi{z)-Gi{-z)^2j2 



oo ^(2fe+l)(o) 



fe=0 



{2k + iy. 



-z 



2k+l 



(8.16) 



matching terms with Equation (|8.15p implies = I and Equation (|8.10p . 

To complete the proof of Proposition l8.4l we now give the proof of Lemma [875l The 
idea is to re-express the polynomial expression for Ri{z) in terms of spherical Hankel 
functions and to achieve simplification by using the Wronskian identity governing h 



(1) 



and h\^\ Recall 



piiz) ^ z'+'h['\z)e 



p\{z) = z' zh^l^'{z) + + iz)h\'\z) 



Substituting into (j8.13p yields 

Ri{z) = 2i \z^+^hf-\z)e-'''{-iy+^z^+^hf-\~z)e'' + z 



(8.17) 
(8.18) 



21 



+ z' {zdh'^piz) + + iz)h^^\z)\ e-^^(-l)'+iz'+i/i|^V2)e' 
+ (-l)'z' \-zdh\'\~z) + + l + iz)hf\-z)] e''z'-+^h^^\z)e- 



2iz' 



{-ir'z^h^^\z)h^^\-z) + l 



+ (-1) 



1 + 121 



z'dh'i^\z)hy>{~z) + {L^- iz)zhY\z)hY\~z) 



,(1) 



+ (-l)'+^z^' z'-dh^i\~z)h\'\z) - {L^+ iz)zhY'{z)hY\~z) 



,(1) 



2iz 



+ (-l)'+iz2' [z^i^dh['\z)h['\~z) + dh['\-z)h['\z)y2u^ 



-2/ 



2i + {-iy+'z' [dh\^\z)h\''{-z) + dhy>{-z)h\''{z) 



= z''[0]=0 . 

(8.19) 
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In the last line we used a Wronskian relation h';^'' {z)dh\'^\z) 



h 



-2iz ^ with the analytic continuation formulae jT] h^'^\z) = (— (— z) and 



dhf\z) 



□ 



This completes the proof of the key Lemma 18.51 and 



thus Proposition [ 

We now complete of the proof of Theorem l6.2l Recall that deformation scattering 
resonance energies satisfy the equation + nG;(eA) = 0. In the (incompressible) 
limit, e — >■ 0, we have the equation }? + riGiifi) = 0, which has two real solutions, 
the Rayleigh frequencies, to which there are associated undamped oscillatory modes 
of the incompressible linearized equations about the spherical bubble equilibrium. 

We are interested in where these Rayleigh frequencies perturb for positive and 
small e. By Theorem 16. 1[ they necessarily perturb to the open lower half plane. 
Propositions 18. II and [ 8.41 enable us to compute, in detail, the location of the Rayleigh 
scattering resonances to which these real frequencies perturb. 



We first consider the case: I = 0. Using properties of Hankel functions (see Ap- 
pendix [DTI} we have Go(eA) = ieX — 1 and the equation for the resonance energies is: 
A^ + zero A — fo = yielding the solution in Equation (|6.2p . This completes the proof 



of Part 1 of Equation (j6.2 



For general ^ > 2, we can, by Propositions 18.31 and 18.41 write the resonance equa- 
tion in the form: 



I 



21+1 



(2^ + 1)!^' 



(2/+1) 



i0)+O(e 



.21+2 



= 0, 



.20) 



A{X) 



where A{X) is even and real-valued for real A and B{X) is odd and purely imaginary 
for real A and r; > is given in (|l.lip . It is clear from equation (|8.20p that the lowest 
order term contributing an imaginary part is at 0{e^^^^). Taking 2Z + 1 derivatives 
in e and setting e = yields 2A(0)A(2'+i) (0) + r,Gf (0)A(0)2'+i = 0. Again, using 
Proposition [831 we find A(2'+i)(0) = -nG(2/+i) (o)A(0)2' and 



QA(2'+i)(0) = [ri{l + l)f {21 + 1)! 



2H\ 



{21)1 



.2^i-^r\+\l + iy{2l + iy. 



n 



{21)1 



.21) 



This proves Part 2 of Theorem [ 

9. Estimation of deformation resonances for I > Kie~^~^, S > 0; 
Proof of Theorem 16.31 . In Theorem l6.1[ we established that there are I + 2 defor- 
mation resonances in the lower half plane. These are solutions of A^ -I- riGi{eX) = 0; 
see Equations (|l.lip and (|6.ip or 

\2 



(6A)^ 



^Wc-\l + 2){l-l) 



hf^ (eA) + eA dh^^ (eA) = 0. 



Equivalently, in terms of standard Hankel functions: 

( -\Hi^\tX) + eX ai/«(eA)=0, 

Ve^We-i (1.2 _ (3/2)2) 2/ " ^ ^ " ^ ' 



(9.1) 



V = 1 + 1/2. (9.2) 
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In this section, we estimate the location of these resonances in the parameter 
regime where the second term dominates. In Section 19.11 we will show that I + 1 
deformation resonances are near the I + 1 zeros of dHl^\eX) by applying Rouchc's 
Theorem. See Appendix IE. 41 for the locations of the zeros of Hankel functions. In 
Section [9.2[ the remaining resonance is found by the method of dominant balance. 

The proof of Theorem 16.31 is based on an application of Rouche's Theorem [2] : 

Suppose f and g are analytic inside and on a simple, closed curve C and satisfy \g\ < 
I/I on C. Then f and f + g have the same number of zeros (counting multiplicities) 
inside C. 



9.1. Proof of Theorem 



Part 1. Let 

vz = eA 



(9.3) 




(i^z)' 



£2 We"' (i/2 _ 9/4) 



i7,(i)(;yz) = 0, V ^1 + 1/2. (9.4) 



We control the zeros of (19.41) in terms of the zeros of dlll^^ 



for / sufficiently large. 
Here, / and g arc defined for the application of Rouchc's Theorem to the locations of 
alH + 2 resonances. 

We prove that the I + 1 roots of Equation (|9.4p . denoted Zj^^i, z^^2 
asymptotically (l large) close to zeros of dH^\vz). This implies, via 
are resonances at the locations: 



z,y^i+i, are 
that there 



A^,(e) 



(9.5) 



Let y'jj^ denote a solution of dH^'{y) = in the fourth quadrant; see Equa- 
tion (|E.17[) . We claim that (/ + g){z) has a unique simple zero, z^^j, in the interior 
of 7^, a small curve about the point z — v^^y^ j. The curve 7^ will be the image of a 
small circle defined in ()9.10p . To prove this, we show that for all z G 7z, 







1 / (^Z)2 






i^z Ve2We"^ (i^2 _ 9/4) 



Hi'\i^z) 



dHl^\vz) 



< 1, 



(9.6) 



for u > Ki{e-^Wef^^^ ^"^^ , where < 9 < 1/3. 

The key estimate, (j9.6p . is a consequence of the following two claims. Below, we 
use C to denote a generic positive constant, independent of z^, e,We. 
Claim 9.1. 

(i) Choose e > and < Ci < C2. Assume ci < \z\ < C2- Then, there exists 
C > 0, depending only on ci and C2, and — 0(1) > such that 



1 

vz 



{vzf 



e2We~' (i/2 - 9/4) 



< 



C 



e2We" 



V > Vq. 



(9.7) 



(ii) There exists > and C3 > such that for all v > and |arg2| < tt with 

kl>C3, 



dHi^\vz) 



Ai{r]) 
dAi(ji) 



where 



(9.8) 
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^ = e2"/3^.2/3^^ ^ / ^ log ^^vi _ yiT^ ; (9.9) 



see also Sections \E.3\ and \E.4\ 

Claim 9.2. Let q G (0, 1/3) and denote by r]'^ the sth solution of dAi{r]) = 0; see 
Section \E. 2[ There exists vi > 0, k > such that if v > vi and s G {1, 2, ...,[(/ + 
1)/2J + 1}, then 



dAi{r]) 



< Ciy^^^+'', on the circle = |?7 : \r] - Vs\ = |s"^/^J^"'} , (9.10) 



chosen to enclose exactly one zero ofdAi; see Equation (|9.22p . 

Inequality (|9.6p is a direct consequence of (|9.7p . (|9.8p . and (|9.10p of Claims [93] 
and ESI 

Assuming Claims IHTTl and l??^ we now prove the assertion of Part 1 of Theorem l6.31 
concerning /+1 of the 1+2 resonances. (The (Z+2)nd is covered by a separate argument 
below; see Section 1^?^ ) 

First, fix the parameters q, vi, and k as in Claim and choose fj G Cn, so that 
the upper bound in (|9.10p holds. For such 77, by the change of variables (|9.9p . z{f}) 
varies in an annulus < ci < \z\ < C2, where < ci < C2 < 00 depend on k and vi. 

By (|9.10p . we write fj G C,, as 

7? = + |s-i/3j/-'?e*'', < 7? < 27r. (9.11) 
Under the mapping 77 1—> Ci (|9.9p . the circle is transformed into the circle C(^: 

C(77) = j.-2/3e-2.,/3 ^ + ^5-i/3^-«g<^) g c^^. (9.12) 

Finally, we map Cq to the curve 7^ in the z-plane using a formula from |22) 

z(C) = 1-2-1/3 C + -^2-2/3 C' + ^C' + C'0(1), C small dEH]) 
Thus, to leading order, 7^ is a still a circle: 

7. = ziCc) = 1 - 2-1/3 Q + O (9.13) 

with radius the same as that of 2~^l'^ Cq, namely 2~^^^ k s~^/^v~'^/^~'' . Since this 
radius is decreasing in both s and v, it can be bounded uniformly. Using this uniform 
bound, we can find fixed ci and C2 defining an annulus which encloses 72 for all values 
of s and v. 

Additionally, if we take y'^, j, to be a zero of dH^ corresponding to 77^ , and approx- 
imated by Equation (jE.17p with k — s, then v~^yls is inside the curve 7^ defined by 
Equation (|9.13p . It is equivalent to show that when mapped into the ( variable, the 
image of v^^y'^, ^ lies in Gc^. From Equation (jE.17p . we see the image of v~^y'^ ^ lies 
in a disk of radius ^ zy-"'/^ which is smaller than the radius of Cc^. With this choice 
of ci and C2, we have, using the definitions in (|9.4p and Claims [9Tl1 and |9^ that 
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Therefore, there exists Ki > such that for any q G (0, 1/3), if ? > /vi (e-^We)^^^^ 
then \ f{z)/ g{z)\ < 1 on the required curve jz- The largest range of Z-values is obtained 

by taking q ^ 5 > and small and thus / > i^i (e^^We)'^^'^^''^ 

Since we have shown dH'i^\vz) ~ has one zero inside ^z, by Rouchc's Theorem 
there is also a single solution, z^^s, of (|9.4p . inside 7^. To finish the proof of Theo- 
rem [531 Pa-rt 1, we compute z^^s and map it to the A-plane. Since z^^s is inside 7^, we 
use Equation (|E.15p evaluated at (|9.12p for the approximation. We find: 

z,., = l+2-i/3g-2W3(;+i/2)-2/3|^^,|^^2-i/3e-'^-'/3^-4/3 |,,;|Vo(5-i/3^-2/3-,y 

(9.14) 

Mapping back to A, using (|9.5p completes the proof. 

Proof of Claim \9.1\ Estimate (|9.7p is straightforward. To prove estimate (j9.8p 
we use the following leading order asymptotics of Hankel functions, uniform in z and 
valid for large v (see Appendix IE.3P 



argz| < TT, 



1 f I- z' 



2\ 1/4 



4C 



dAi{T^) 

7/2/3 ■ 



argz < TT. 



In particular, for any fixed < (5 < 1 and sufficiently large u 



Hi^\vz) 



dH^}\vz) 



< 



1-S 



.1/3 



1-zJ 



1/2 



Ai{ri) 



dAi{-q) 



(9.15) 



To show this is bounded near z 1, we use the Taylor series of Q{z) about z = 1, as 
given in Appendix IE. 31 



so, due to the boundedness of z. 



C = -e"/32i/3(z - 1) [1 + O (z - 1)] , 
1/2 



(9.16) 



C 



< C. This completes the proof of 



and therewith Claim I^TTI □ 

Proof of Claim Since dAiijj'^) = 0, Taylor expansion of Ai and dAi about 
rj = rf^ implies the existence of analytic functions Ai and A2 , such that 



A^ir,)=A^i,i)+A,i7^)i7^-,if 
dMir,) = v'sMri'sKv - v's) + Mv){v - v'sf- 

To obtain (|9.18p we also used that d^Ai(r]) = rjAi{r]). It follows that 

Ai{ri) Ai«)+Ai(?7)(?7- 7/^)2 _ 1 



dAi{-q) r/^Ai{rf^){r] - T]'^) + A2{r]){r] - i]'^)"^ v'^ir] - v's) 



(9.17) 
(9.18) 



(9.19) 



Thus, we now estimate a'^(r]). Using the domain specified in Equation (|9.10p . we have 
the bound 



1 



(9.20) 
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Wc complete the proof of Claim by showing |a'j,(?7)| < C on a circle around if^ that 
is small enough such that it encloses no other dAi{j]) zero. The specific circle used is 
defined by 

177 - 77^1 = Cs-i/3t.-9 

and is determined by the spacing between consecutive zeros of dAi(ri); sec Equa- 
tions (|9.22p and (|9.27p . There are two cases to study: 

• First, consider the case s > s'l > s'q > 1 (sg and s[ independent of I). 
Using (jE.8[) . giving the locations of the larger dAi{r]) zeros, {rj'J\ s > Sq}, we estimate 
the spacing between the zeros as 

Ws^,-v's\ = l{'-Tf'^-'/' + 0{s-'/') (9.21) 
so that the spacing satisfies the lower bound 

Ws+1-V's\>^s-'^', ^>So (9-22) 

where k = (37r/2)2/3/3 w 0.93693. So, first, we consider the open disk \i] — rfj < 
^s~^^^, which clearly contains only one zero of dAi. 

Now we estimate the Ai and A2, appearing in (|9.19p . Using the Cauchy integral 
formula and Ai (77) obtained from (|9.17p we have 

(9.23) 



Using the asymptotic expansion of Ai{ri) on the negative real axis from Appendix lE.il 
and the location of the t]'^ from (jE.7[) we have 

Cqs-^/^ < sup \Ai{ij)\ < ClS-l/^ s > s[. (9.24) 



This allows the estimate for I?? — ?7^| < ks' 



-1/3. 



, . / M 1 2\Ai{i^)\ r 1 , , , 
|Ai(r;)|<— sup f :\dt\, 

S-l/6 f 1 

-'^1 m^f r, rl^^l' by equation dSSl, 

(ks ^I-^Y J|t_,,^|=Ks-i/3 \t-'q\ 

s^i/^ 1 / 
- ^1 — ^1/3^2 — ^173 \ Ti r 1 1^*1' by the triangle inequality, 

(,KS ' j KS \rj ?7g| J|t-r;^|=/ts-l/3 



1/3 1/3 _ _ I 

Next, we apply the last inequality to the smaller circle \r] — ?;^| = Cs^^^^v^'^ we can 
estimate 

\Ai{v)\ < C — j < Cs'/^ ^ _ < Csi/2. (9.25) 

[KS -^/■^j^ 1 — \ri — rig\ / [ks '-'■^) 1 — Lv 1 
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Similarly, 

\A2{v)\ < Cs'/^ 

The full numerator of (|9.19l) can be estimated, using (I9.24[) and (|9.25p . 



1 + A^{^i)-^A,mr1 < 1 + IMi^'Jf' \A,{r^)\ \v - v'sf 

and 



For the denominator, using (|E 



< 



'7s 



-^\A^{^,)\-^\AM\\^l-r^. 



(9.26) 



< 1 + Cz/ 



-2q 



l + [i,Ai{i^)] 'A^ivKv-rfs 



< 1 + Cv Combining these esti- 



Thereforc, 

mates, we find \a'^,\ < C. 

• The case 1 < s < s'^ is handled similarly. Since we do not have an explicit 
approximate formula for the zeros of dAi zeros, we introduce the among the first s'l 
zeros 



26'^ 



1<S<S, 



'7s 



For 1 < s < s[, define Af and A| via ((9Tfl) and (pl8| . Then, 

, _ l + ^z(77^)-iAf(r;)(r7-r7;)^ 
l + r;^^z(770-iA^(77)(77 - 770' 



(9.27) 



(9.28) 



In the disk I?? — 77^] < 6'i/~'^ (in the theorem Cs~^/^ = S') wc use the Cauchy integral 
formula to make the estimates |Af (t;)! < C and |^2('/)l — ^7 which give |a^| < C. 
This completes the proof of Claim 19.21 and therewith. Part 1 of Theorem 16.31 □ 



9.2. Proof of Theorem [6T3] Part 2. We now discuss the single resonance, away 
from the arc. As before, we begin by seeking solutions of Equation ()6.1|) for large I 
using the asymptotics in Equations (jE.9P and (jE.lOj) . However, we now seek zeros in a 
domain not associated with zeros ofdAi{ri). In particular, we focus where \r]\ is large 
and largjyl < tt. Using the leading term of the asymptotic approximation to Airy 
functions in Equations (|E.1[) and (|E.2[) and Equation (|9.9p we obtain the approximate 
equation 



2 2 
v z 



2 We 



We 



1 



2\l/2 



(9.29) 



where A = i^z/e. For v 3> 3/2, we can drop the terms involving 9/4, square both 
sides, and write an approximate polynomial equation for the resonances in z: 



4, 24/'1^2 2 2 / 1 \ 2 24/'l^ 



0. 



(9.30) 



All neglected terms are of higher order in e and will be absorbed into the final error. 

o 2 

set a = and obtain: 



We identify the large parameter 

z^ + N^z^ - 



€ f 

Wc ' 



2 2 

a z 



Wo 

0. 



(9.31) 
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In the regime, we have N > K2, and the dominant balance is between 

the first two terms. Choosing a solution in the lower half plane leads to z w —iN. 
Accordingly, we look for a solution in the form 



z{N) = -iN + Z0 + 2-i7V-i + z-2N-'^ + 0{N-^) 



(9.32) 



Iteratively, we compute zq — 0, z_i ~ i{a — l)/2, and z_2 = 0. Changing variables 
back to A leads to the desired result. This completes the proof of Theorem [ 

Appendix A. Spherical coordinates. 
Spherical coordinates in M.^: 



X = r sm a cos ( 
Standard orthonormal frame 

e = 





y = r sm sm i 



-der, cf> 



rcos9 




grad, div, A in spherical coordinates 



1 df 



grad / = Vf=^r + -%e , . 

or r oO r sm t/ 0(f 

divF = V.F=l|-(r2F,,) + ^^ 
or r sm a at) 

1 







(sin^^^e) 



1 



r sin 6 1 



-Fa 



div grad/ = A/ = A,/ + ^Ag/ 



A. / 
As / 



1 

sin^^ 



a/ 

9r 



radial Laplacian 



sin 9 



dl 

de 



1 

sin^e 902 



Laplacian on S 



Surface and Volume elements: dCl = sm9 dO i 



dx = r sin dr dO i 



(A.l) 
d^i (A.2) 

(A.3) 
(A.4) 

(A.5) 
(A.6) 



Appendix B. Spherical Harmonics. y,"(6',(^) = y,™(0), Z > 0, |m| < ^ 
denote the spherical harmonics. These are eigenfunctions of the Laplacian on S^: 



AsV(^,, 



l{l + l)Yr{e.cj,). 



(B.l) 



The functions {i^/"} form a complete orthonormal system in L^{S^): 
(Yj"*, Y;i"')i2(52) = (5/;',mm' with L^{S^) inner product: 



27r /'IT 



(./i.9)l2(s2) — 
We shall make use of the following explicit formulae 



f{9, (/)) g{0, (/)) sin 9 dd dcj) ^ j f{n) g{n) dfl 



J- n 



(B.2) 

(B.3) 
(B.4) 
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Appendix C. PDEs for a bubble in a compressible fluid and lineariza- 
tion. The inviscid, compressible fluid phase is governed by the Euler and conservation 
of mass equations with normal velocity and stress boundary conditions: 

= dtu+ {u-V)u+-Vp(p), xeM.^\B(t) (C.la) 
P 

= dtp + V -{pu), xeM.^\B{t) (C.lb) 

u-n = 5tR-n, dB{t) (C.lc) 

PbubbiclaB(t) - -Pfluidlas(i) = • dB{t). (C.ld) 

The pressure within the fluid is assumed to obey an equation of state: p = pnuid = 
p{p). The pressure of the gas within the bubble is assumed to be constant and to vary 
inversely with bubble-volume: 

k 

Pbubbic = -Pb(R.) = (^-2) 

For air (composed of mostly diatomic gases), the ratio of specific heats 7 = 1.4. Equa- 
tion (|C.ldp is called the Young-Laplace equation [15] and states that the discontinuity 
in pressure at the interface is proportional to the mean curvature H due to the rela- 
tion V • fi = 2H; see Lamb [13] (Article 275, Equation 5), or do Carmo [7] (Section 
3-3, Example 5). 

Equations (jC.l[) have a spherically symmetric equilibrium solution: 

u = 0, Il{e,(t)) = ar{e,(j)) (C.3) 
p = Poo, P = Poo, Pb = Pcq (C.4) 

where the equilibrium bubble radius, a and pressure, Pcq are given by (jl.Sp . 

To make the fluid compressibility more explicit, we deflne a variable sound speed 
= ^ and rewrite the continuity equation ()C.lb[) in terms of the pressure, p = p{p): 

dtp+'Vp-u + c^pV • u = 0. (C.5) 

C.l. Nondimensionalization of the compressible Euler equations. We 

shall use the equilibrium bubble radius, a, as a length scale, and a typical bubble 
wall velocity, U, as a velocity scale to nondimensionalize the compressible Euler equa- 
tions (jC.ip . In this section, we denote the nondimensional variables with subscript 
asterisks. Thus, the non-dimensional spatial coordinate is denoted x.^. Now, introduce 
non-dimensional variables: 

B = aB^ x = ax* t=^t^ u = U p = PooU'^ P* (C6) 

dB = adB, Vx = iVx. dt = -dt, p = Poo P* Pb = P^U^ P*b, (C.7) 
a a 

In terms of these non-dimensional variables, the first two equations of the Euler system 
become: 

9t,u« (u* • Vxju, + — Vx.p*(p*) = 0, x^eM.^\B^{U) (C.8) 
p* 

^t.p* + Vx.- (/9*u,) =0, x,eR^\B,{U). (C.9) 
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Remark C.l. // Vx, A u^(t = 0,x*) = 0, i.e. the initial velocity field is 
irrotational, it follows that Vx, A u.^,{t,x.^,) = /or all t; see JIS}/ . In this case 
UH,(t,x*) = Vx, <^*(t,x*). Alternatively, we may write the continuity equation (|C.9[) 
in terms of the non-dimensional pressure, p*, as: 

dt,P* + Vx,p* • + c* V* Vx, • u, = 0. (C.IO) 

Here, we have introduced the non-dimensional speed: = iff)'^ — ^dp^iP*)- Note 
that , u* , and p* are all functions of (x* , ) . 

Next, we nondimensionalize the boundary conditions. Using the substitutions 
above we have 

u* • n = (?t,R* • fi, x*e9i3*(t*), (Cll) 

ns. -P* = Vx,-n, x,edB,{t,) (C.12) 

We 

where We denotes the Weber number, denoted by 

^''^ a ~ CURVATURE- ^^'^"^^ 

We now express the dynamic boundary condition more explicitly. First note \B{t)\ = 
\aB^,{t^,)\ = a^\B^,{t^,)\. Furthermore, the bubble pressure at equilibrium is given by 

Peq = k^Peg (fo^)^ =^ PB(t) ^ Peq (^)^ ThuS, USiug 

the expression for P^q in (|1.5p . we have 



/OooC/2 j V|5*(i*)l/ V 2 Wcy 

Here, Ca denotes the Cavitation number 

_ Poo _ EXTERNAL PRESSURE 

~ ~ KINETIC ENERGY PER VOLUME' ^ ' ' 

The nondimensional system is given by Equations (jC.8|) . (jC.9[) . (jC.ll|) and (jC.12[) . 

Finally, we conclude this section by displaying the spherical bubble equilibrium 
solution in non-dimensional variables: 

u, =0, R, - v[e, 0), p, ica, - 1, P,B, = ^ + ^ (C-16) 

C.2. Nondimensional Euler equations in Center of Mass Coordinates. 

The bubble's dimensional center of mass is defined 



«_(t):= / xdx (C.17) 

JB(t) 

and is nondimensionalized by ^cm(0 ~ '^^cm *('•*)■ this section, variables and 
operators defined with respect to the moving coordinate system are denoted with a 
prime ('). Let the moving coordinates be defined 

x:(t,):=x,-|,„,(i,) (C.18) 



40 



A. M. SHAPIRO AND M. I. WEINSTEIN 



such that 

dx; = 0. (C.19) 



's:(t.) 

We will use the following substitutions. 

Vx. = , n = n', R, = + (t, ) 
u4x*,t*) = u'Jx* - $^^^{t^),t^) = u'Jx'^,t*) 

p4^,,Q = p:i^:,Q, (C.20) 

(x, p', (xt ,U)P, B. (R* .t,) = PlB', {K , i* ) , (x, , ) = 4 (xt , ) . 

The resulting system of nonlinear equations is: 

= {dt,-dtA.m* ■ VxiX + K- Vx'J K + ^Vx'K, x^eM^BU^,) (C.21a) 

^ = {dt,-dtAc^..-^K)p'*+"^KP'*-<+Cp'* Vx'.-<, x:eM3\B:(i,) (C.21b) 
u:-n' = 5t.(R:+^,^,)-n', e dB^U) (C.21c) 

n's: -P'* = 7^ Vx; • n', x: G 9B:(t,) (C.21d) 

We 

= / x', dx',. (C.21e) 

C.3. Linearization of Euler equations about the spherical equilibrium. 

To simplify the notation, we drop asterisks and primes in Equations (jC.2ip . and 
linearize about the nondimensional equilibrium (|C.16|) by taking: 

u = + 5ui+e'(52), R(6i,0,t) = i?(6i, 0, t) f = [l + (5i?i(6i,0,t)+e'(52)]f 

(C.22) 

Ca 

p= — + 5pi + 0{5^)- p = pip), p^l + 5pi + 0{S^) (C.23) 

lcm(i) = + HcmAt) + 0{S') , = a,.(Ox + ^l,ym + 6,. (C.24) 



Consider initial conditions that are perturbations from equilibrium of size 6, 5 < I. 
We expand quantities in Equations (|C.21[) in powers of S: = ^(1 + ^Pi + • • ■ ) = 
+ i5 ci^ + 0[S^), where M denotes the Much number 



M= ('^(1)^ ' = — 

\dp J Ceo 



INERTIA 



COMPRESSIBILITY 



1/2 



(C.25) 



In addition to the above expansions in powers of 6, we require the implied expan- 
sions for h and div n. From the relation 



F{r,0,(j),t) = r - R{e,(j),t) = 
which defines the bubble surface, dB{t), we have using the polar coordinate represen- 
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tations, (|A.3p and (|A.4[) . for grad and div 



n 



r=R 



grad F 
I grad F\ 



i-Tl9eR0-nked,R^ 



r=R 



fl2 



{dgRy 



fl2 sin^ 



(C.26) 



1-r-S [ deRiB + -±^d^Ri^ ) + 0{S^) , 



r=R (C.27) 



(div n)|^^^ ^ [div i^l-r - S deRiO - S ^ + 0{S^ 

= 2-6 i2 + As)Ri+0{6^) 
It follows that 
u • n = +S drUi + 0{S^) 

sin0cos(/) + 9t^i,y (t) sin sin (/) + dt£,i.z{t) cosO ) + 0{S ) 
9tR-n = +SdtRi + 0{S^) 

For the bubble pressure, we use the following expansion of the volume, 



B{t)\ = del) dO sine 

Jo Jo J 



r dr = —TT + 5 



Ri{e,(j),t) sin d de + o{s^) 



Therefore, 

471- \ 7 



\B{t)\ 



= 1 - 37 <5 {Y,',R,)l^s-) + 0{S^), where ¥,"{9, . 



1 /I 

2 V tt' 



The center of mass constraint can be expanded is follows: 



xdx 



Bit) 



de 



R{e,4>,t) 



r r{9, 0) sinfl dr 

- / R'^{9,(P,t) sin6lsin0 sine* dd d(j) 

Jo Jo \ cos9 I 



(C.28) 



' sm t/ cos ( 



Since R = 1 + S Ri +0(5^), we have at order 6: 

j-TT /sin COS , 

/ / i?i(6',0,i) sin^sin^ sin6'd6'd0 = O 
Jo Jo \ ^^ofl / 



(C.29) 



COSf 



Equivalently, we can express these three orthogonality conditions in terms of spher- 
ical harmonics of degree I = 1: Yi^0,(l)), Yi°(6l,0) and Yf\0,<l)). Indeed, by Equa- 
tion (|B.4p 

n6'cos0 = yf^-yi\ y^sine'sint/) = i (yf ^ +r/) , ^^J^cos9 = Y^. 

(C.30) 
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Therefore, the three orthogonahty conditions of Equation (jC.29P arc equivalent to 
Yi™(-))2^2(52) ~ 0, m ~ —1,0,1. We can also re-express 9t^j.m ' given 
above, in terms of as follows 

= ^5 {V2dtS,iAi) {Y\-Yl)+V2idti^,y{t) (Ff 1 +F/)+2a^el,.(^)n")+^^('5')• 
(C.31) 

Substitution of these expansions into Equations (jC.21[) yields the following equations 
for the linearized perturbations to the fluid phase: 

= 9tUi + Vpi, {) = dtPi + ^V-Mi (C.32) 

Assuming the flow is irrotational, we have u = 5 Ui + 5^ U2 + 0{5'^) = 
6 Vtfi + (5^ Vip2 + 0{6^). Since ui = Vtpi, we have = V{dtipi+pi) and 
therefore 

= dt(pi+pi, O^dtPi + ^A^i 
implying the wave equation for tpi: 

M^a^Vi - =0, r > 1, t > (C.33) 
The linearized boundary conditions and center of mass constraint are 

drifi = dtRi + (9tlem ' f) {9, <t>,t), r = 1, t > (C.34) 

dt^i^ il {^ + ^^ {YS,Ri)Y^^^^{2 + As)Ri, r = l, t>0 (C.35) 
=0, |m|<l (C.36) 

As initial conditions we take: 

i?i(6', 0, i)|j^Q sufficiently smooth 

(pi(r,0,0,Olt=o' 5t¥'i(r-,0,</',t)lt=o =0, r->l, (C.37) 

Proposition C.l. Let ipi{r,0,<^,t), 4cmi(*) denote a sufficiently 

regular solution of the initial-boundary value problem for the wave equation (|C.33P on 
the region ?■ > 1 with linearized kinematic and dynamic boundary conditions (jC.34p 
and (jC.35p . linearized coordinate normalization (center of mass) condition, and initial 
conditions (jC.37p . Then, 

^e™,iW=0, t>0. (C.38) 

It follows that the linearized perturbation: = ipi, (3 = Ri satisfies the initial- 
boundary value problem (jl.Sp . (|1.9p of the Introduction. 

To prove Proposition IC-H we first introduce the projection operators onto the 
|m| = 1 spherical harmonics 

vr = {-,Yr)ms^)Yr. (c.39) 
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Since Vt^ commutes with dt and A, for any \m\ < 1, the function U — VT^dt^pi 
satisfies the wave equation 

{M^d'i - A) C/ 0, r > f (C.40) 

Moreover, applying to the dynamic boundary condition (jC.35p and using (|C.36p 
yields the Dirichlct boundary condition 

C/ = 0, r = f. (C.4f) 

ft follows that U = VY^dtipi is identically zero on r > 1. Therefore, 

Since Vi^ipi is a time-independent solution of the wave equation for r > 1, clearly 
AV^ifi = 0, and therefore (A^ - 2r-^)F{r) = 0. Thus, F{r) = ci ii(r) + ca yiir) for 
some constants Ci , C2 and therefore 

rrVi{r,0,<l>,t) = ic,nir) + C2yi{r)) rne,^), and 

Evaluating (jC.42[) at r = 1 and t ^ we have, using the initial condition ipi{r, 0,(j),t = 
0) = for r > 1, 

cu'i(l) +C2yi(l) =0, ci(9^ji(l) + C2a,.yi(l) =0. 

implying ci = C2 = 0, by linear independence of {ji{r),yi{r)}. Now applying T'™ to 
the kinematic boundary condition (|C.34[) and using (|C.36[) . we obtain 

( W • f) = 0, \m\ < 1 (C.43) 

We now claim that <9t^cm.i(*) — 0, ^ > and therefore $cm.i(*) — 0, t > 0. To 
see this, observe that 

t^t^cm,i(0-r = 9t^i:E(i)sin6'cos6' + (9t^iy(t) sin6'sin9!i + dt^iz{t)cos9 

Taking inner product with Y"\ \m\ < 1, yields dt£,ix{t) = dt£,iy{t) — dt£,iz{t) = 0, 
t > 0. This completes the proof. 

Appendix D. Hankel Functions. Separation of variables, applied to the to the 
wave equation in M.^, in spherical coordinates, leads to the spherical Bessel's equation 

r^R"{r) + 2rR'{r) + (r^ - 1(1 + 1)) R{r) = 0. (D.l) 

Two linearly independent solution are the spherical Bessel functions 

Mr) = y^J/+i/2(r), yiir) = ^Yi+./^ir) (D.2) 

where Ji+1/2 and Y;+i/2 are Bessel functions of order I + 1/2. The outgoing spherical 
Hankel function h^^\z), z e C, is obtained by taking a linear combination of ji and 

yi- 

h'i^\z)=ji{z) + iyi{z). (D.3) 
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In terms of the Hankel functions, H^^^^2{z) = Ji+i/2{z) + i^;+i/2(^) and 

D.l. Polynomial Representation. We can represent the spherical Hankel 
functions as in [JD] by 

h']'\z)^z-'-^pi{z)e' (D.5) 
where pi{z) is the polynomial of order / given by 

k=0 ^ ^ ^ ' n=0 ^ ' ri=0 

D.2. Analytic continuation. The Hankel functions have the following analytic 
continuations and symmetries 



(ze-) = -e— ^(1) (-) ^ ^(2) 
(ze-"*) = -e-^^'i/^i) (z) ^(2) (-^ ^ ^(1) ^ ^e^^'Hl^^ (ze") (D.8) 

D.3. Derivatives. The derivative of a spherical Hankel function can be ex- 
pressed in terms of spherical Hankel functions in a variety of ways: 

These relations are also valid for ji , y; , /i^ . Similarly, 

^ - = ^ (D.IO) 

These relations are also valid for the Bessel functions Ji, Yi, and the incoming Hankel 



(2) 



function 

Finally, a limit we require is a consequence of the above recursions and the poly- 
nomial representation of /i^^^: 

Appendix E. Asymptotic Expansions of Special Functions. 

E.l. Airy function asymptotics. The Airy function and its derivative can be 
approximated for large argument by (Olver |23| ) 

Ai(T) = i^-i/V-i/^e-M + ^ |argr|<7r (E.l) 

dAt{T) = -^TT-'^\^^^e-^[l + 0{\fi\-^)], |argr|<7r (E.2) 

Ai{-T) =7r-i/V~i/''cos(/x-7r/4) [l -I- tan(/x - 7r/4)0(/i-i)] , |argr| < 2tt/3 

(E.3) 

Ai'{-T) = Tx-^I'^T^I'^ cos(/i - 37r/4) [l + tan(/i - 37r/4)0(/i"^)] , |argr| < 27r/3 

(E.4) 
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where /i = ^t^^^. Additionally, we have the inequalities 



|^z(r)|<c(l + MV4)"' 
|c>^z(t)| <c(l + |r|i/4) 



e 

-Ml 



|argr| < tt 
larErl < tt. 



(E.5) 
(E.6) 



E.2. Zeros of Airy Functions. All of the zeros of the Airy function and its 
derivative lie on the negative real axis. The larger zeros have the following approxi- 
mations [nmjini. 

Asymptotics for zeros of Ai{z): 



-7r(4s- 1) 



2/3 



So, ■ 



/ + 1 



Asymptotics for zeros ofdAi{z): 



Vs 



-7r(4s - 3) 



2/3 



'0, ■ 



(E.7) 



(E.; 



E.3. Hankel function asymptotics. The Hankel functions can be approxi- 
mated for large order v, uniformly in z in the sector |arg2;| < tt by 



ai/(^)M^-4e-/3±(^^j ^' (E.IO) 
where 

V= ^r^,i^ ,„gi±^ - (E.n) 

6 t z 

We note the useful Taylor series 

^C'/2 = (-log(z) + log(2)-l) + ^ + 0(z4) asz->0, (E.12) 

= -^*V2(z-lf/2-t-Mi^il!^! + 0((z- 1)7/2) asz^l, (E.13) 



5V2 

iiT i i 1 / 1 
y 2z 24z 



= -^z+------+0(-] asz^oo. (E.14) 



which yield the limits lim2__j.o C(^) = 00, lim^_j.i C(z) = 0, lim2_j.oo Ci^) — oo- For ( 
close to zero. IE .iT] is inverted to yield 

z(C) = 1 - 2-1/3 ^ + ^ + C'Oil). (E.15) 

E.4. Zeros of Hankel functions. Concerning the zeros of Hankel functions a 
consequence of the analysis in [35] is the following: 

Theorem E.l. Let v = I + \/2, where I > 2 is an integer. The I + 1 zeros of 
dHi^\y) and the I zeros of Hi^\y) lie near arcs in the lower half plane . 
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The detailed locations of the previous zeros given in the following 
Proposition E.2. The I zeros of Hl^\y) are 



(E.16) 

for k = 1, . . . , \_{l + 1)/2J , where rjk is the k-th zero of Ai. 
The I + 1 zeros of dHi^\y) are 

vIm = - + 2-l/3e-2../3^-2/3 |,^, I + =|.2-2/3^-4W3,-4/3 |2 + 0(,-4/3y 

(E.17) 

for fc = 1, . . . , [(Z + 1)/2J + 1 where rj'f, is the k-th zero of dAi. 

We comment on approximating the zeros oidHl^\vz) by Equation (|ETf)) . which 
is closely related to the analysis of Olver [23 • The proof is analogous for hI^\uz). 
Introduce the variables 

C = C{z), ri = 77(C) = e^'^'/^i/^/^C; see Equation M7n\ . (E.18) 



Using the uniform asymptotic approximation for dHy [vz) in (|E.10p we can 
rewrite the equation dH^\y) = 0, for = / + 1/2 large as 

dAi{rj) + £-^Az(r7)g [1 + 0{u-^)] = 0. (E.19) 

For V large, the zeros, 77, are well approximated by zeros, {77^}, of 9^7(77); see (jE.Sp . 
These lie on the negative real axis; see Figure IEJI (left). Then, the asymptotics of 




Fig. E.l. Complex mappings of some dAi{ri) zeros (left plot) to approximate zeros of 

dHl^\vz) (right plot) by tra nsformations 77 1— >■ i— >■ 2 vz. The first map, 77 1— ^ isgiven 
by inverting Equation SE.lSd . The second map, t-^ z, is computed using Equation SE.lSd . 



Airy functions and the properties of the mapping rj t-^ 1-^ z 1-^ vz can be used to 
show give approximations to the zeros of dH^{y) and to establish their location near 
an arc in the lower half plane. 

Appendix F. Proof of Proposition 15.21 In this section we prove Proposi- 
tion [5?2l which states: 

Res -2 = OL-'/^) . (F.l) 

A=A,,, A2 +riGi{eX) \ J 

The main tool of this analysis is the uniform asymptotics of Hankcl functions. 
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Proof. To compute the residues, in (|5.3p . first caneel a factor of A, use the differ- 
ential equation satisfied by spherical Hankel functions and relations in Appendix |D] 
to obtain: 



Res — '^^ , ^, = Res —, = lim 7 

A=A,,, A2 + nG, eA) \=\i., , , _ h^' (eA) A^A,, / ,,(1)' 



- lim 

A^A, 



1 + -1 



(6A)2 



eA 



-1 i?i'''(eA) 



ff(i)(eA) 



(eA) 

-1 7?i'''(eA) 



2eA ij(i)(eA) 



(F.2) 



In particular, we have used 



v^l + 1/2. 



Now we expand the terms in the last line of (jF.2[) about A;.j for large I. We recall 
the approximation for these resonances from Theorem 16.31 



1 + 1/2 



1 - 2-i/3(? + l/2)-2/3e-2W3,^^ + 



, s = l,2,.. 



(F.3) 

There will be two cases to consider: resonances corresponding to Airy prime zeros 
•q'g for s = 1, . . . ,so which are 0(1) and for s w 1/2 which we will see are 0{P/'^). 
Specifically, 



Btt 



2/3 



(4s - 3)2/3 + 



S= So,... 



l + l 



+ 1 



s w 1/2 — > 00. 



Since 



we have 



{y-l/2){u+l/2) 



1(1 + 1) _ 

(eA)2 ^ 1/2 [1 _ e-2^V32-i/3i/-2/3^^] 

Z(/ + l) [1 + 0(1^-2/3), s small, 



(eA)2 0(1), 



s large. 



(F.4) 



(F.5) 



(F.6) 



For the Hankel terms in (jF.2p , letting eA ~ vz and using Appendix IE.3| we have 

g^'^V^) ^ e2W3i ('i^^ ^,-1/3 ^ ^ ,2W3,2/3. 

i?i'^(i.z) Ac; ^*(t) 

We simplify this ratio for 2 « 1 using the Taylor expansion (jE.13[) . which yields 

C = 2i/3e-"/3(z - 1) + 0((z - 1)2) . (F.7) 
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Thus, C"i(l - z2) = _(i + 



l + z 



22/3g-2../3 + . . . and 



2l/3e-x>/3 

g27rt/3j^2/3(A _ g27ri/3y^^ + ' ' ' , where the dots indicate higher order terms. Now for 
the Airy functions ^^^^^J^ ~ ^(1) for smaU s. From Appendix IE. II 



dAi{T) 

Ai{T) 



T large, |argr| < tt. 



So, for large s 

fdMj)\ 

Combining results 



( dAijT) 
\Mt) 



-TTi/3 



3tt 



2/3 



2/3 



1/2 



-■Ki/6 



1/3 



(F.; 



,1/3 



Oil), 



s small, 



,=,2W3,2/3^ l^e- V6 (f 1.1/3 + 0(i.-2/3)^ , ^ 1/2 large. 
Wc can now conclude the proof. For s large 



(F.9) 



2eA Hi'\eX) 



Thus 
Res 



lim 



A=Ai,3 A2 + riGi{eX) a^a,,. 



1 + t'n -I + - -0{l) - 



(F.IO) 



(F.ll) 



For s small, 
A 



Res 



= lim 



(F.12) 



A=A,., A2 + riGi{e\) a^a,,. 
The asymptotic formula (|F.ip is a consequence of (|F.11[) and (|F.12[) . □ 
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